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Perturbations of black holes, initially considered in the context of possible observations of astro- 
physical effects, have been studied for the past ten years in string theory, brane-world models 
and quantum gravity. Through the famous gauge/gravity duality, proper oscillations of perturbed 
black holes, called quasinormal modes (QNMs), allow for the description of the hydrodynamic 
regime in the dual finite temperature field theory at strong coupling, which can be used to pre- 
dict the behavior of quark-gluon plasmas in the nonperturbative regime. On the other hand, the 
brane-world scenarios assume the existence of extra dimensions in nature, so that multidimen- 
sional black holes can be formed in a laboratory experiment. All this stimulated active research 
in the field of perturbations of higher-dimensional black holes and branes during recent years. In 
this review recent achievements on various aspects of black hole perturbations are discussed such 
as decoupling of variables in the perturbation equations, quasinormal modes (with special em- 
phasis on various numerical and analytical methods of calculations), late-time tails, gravitational 
stability, AdS/CFT interpretation of quasinormal modes, and holographic superconductors. We 
also touch on state-of-the-art observational possibilities for detecting quasinormal modes of black 
holes. 
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I. INTRODUCTION: NEW ASPECTS OF BLACK HOLE 
PERTURBATIONS 

Isolated black holes in equilibrium are intrinsically sim- 
ple objects: They are much simpler than ordinary matter 
around us because they can be described by only a few 
parameters such as their mass, angular momentum, and 
charge. This simplicity, however, is apparent only be- 
cause one never has an isolated black hole. Black holes 
in centers of galaxies or intermediate-mass black holes al- 
ways have complex distributions of matter around them, 
such as galactic nuclei, accretion disks, strong magnetic 
fields, other stars, planets, etc., and are therefore ac- 
tively interacting with their surroundings. However, even 
having been removed from all macroscopic objects and 
fields in space, a black hole will interact with the vacuum 
around it, creating pairs of particles and evaporating due 
to Hawking radiation. 

Thus, a real black hole can never be fully described by 
its basic parameters and is always in the perturbed state. 
Whenever you want to know something about stability, 
gravitational waves, Hawking evaporation of black holes, 
or about the interaction of black holes with their astro- 
physical environment, you have to start from an analysis 
of their perturbations. Moreover, the first moments after 
the formation of a black hole due to gravitational col- 
lapse of matter, the black hole is in a perturbed state. A 
perturbed black hole can be described by the metric g^ 
that can be written as 

9nv = 9% + 59w (LI) 

The metric g^ is the space-time of the nonperturbed 
black hole when all perturbations have been damped. It 
is called the background. In the lowest, linear approxi- 
mation, the perturbations 8g^ u are supposed to be much 
less than the background Sg^ <C g® iu . For example, g^ 
can be the well-known Schwarzschild or Kerr solutions. 



Once a black hole is perturbed it responds to pertur- 
bations by emitting gravitational waves whose evolution 
in time can be conditionally divided into three stages: 

1) a relatively short period of initial outburst of radi- 
ation, 

2) usually a long period of damping proper oscil- 
lations, dominated by the so-called quasinormal 
modes ( "quasi" here means that the system is open 
and loses energy through gravitational radiation), 

3) at very large time the quasinormal modes are sup- 
pressed by power-law or exponential late-time tails. 

The major part of this review is devoted to the sec- 
ond stage of the evolution of perturbations represented 
by quasinormal modes. There are a number of reasons 
for such special interest in quasinormal modes. First, if a 
new generation of gravitational antennas such as LIGO, 
VIRGO, LISA, and others manage to detect a gravita- 
tional signal from black holes, the dominating contribu- 
tion to such a signal will be the quasinormal mode with 
lowest frequency: the fundamental mode. 

Another motivation comes from the famous dual- 
ity between supergravity in anti-de Sitter space-time 
(AdS) and conformal field theory (CFT), AdS/CFT 
correspondence^. The AdS/CFT correspondence says 
that quasinormal modes of a (D+l)-dimensional asymp- 
totically AdS black hole or brane are poles of the retarded 
Green's function in the dual conformal field theory in D 
dimensions at strong coupling. The latter allows us to de- 
scribe various properties of strongly coupled quark-gluon 
plasmas which cannot be studied by traditional pertur- 
bative methods of quantum field theory. Recent achieve- 
ment in this direction was the conjecture of the universal 
value 1/A-k for the ratio of viscosity to the entropy den- 
sity in quark-gluon plasma 3 . The universal value 1/4tt, 
which follows from string theory, is the same for a variety 
of dual gravitational backgrounds and is close to the one 
observed at the Relativistic Heavy Ion Collider (RHIC) 4 . 

Usually we can talk about the following types of black 
holes: 

1) supermassive black holes in centers of galaxies with 
masses M ~ (10 5 - 1O 9 )M , 

2) intermediate mass black holes M ~ 10 3 M Q , 

3) stellar mass black holes M ~ 10M©, 

4) miniature black holes (intrinsically higher dimen- 
sional in string theory and brane-world scenarios) 
with masses of at least a few Planck masses or more. 

Perturbations of the first three types of black holes 
are important for astrophysical observations, while the 
fourth type gives its own motivation for study. It is 
worth remembering that in D — 4 space-time, due to the 
uniqueness theorem, there is only one solution of the Ein- 
stein equations that describes a stationary rotating black 
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hole - the Kerr solution. Thus, we have only one can- 
didate for a black hole metric in four dimensions, which 
is the Kerr metric. In D > 4 space-times the situation 
is different because there is no usual uniqueness theorem 
and, as a consequence, there are a number of "black" 
solutions with an event horizon such as black strings, 
branes, rings, saturns, etc£. These solutions have differ- 
ent topologies, for instance a black ring has the topology 
of a doughnut, and many of them are probably grav- 
itationally unstable. In order to learn which of these 
higher-dimensional solutions are stable, and thus could 
exist in nature, one needs to find their spectra of gravita- 
tional quasinormal modes. Stability is guaranteed when 
all quasinormal modes are damped. Stability of higher- 
dimensional AdS black holes is also important because 
the onset of instability of a black hole corresponds to the 
thermodynamic phase transition in the dual field theory^. 
Thus, quasinormal modes of higher-dimensional black 
holes are connected with stability, thermodynamic phase 
transitions and the hydrodynamic regime of strongly cou- 
pled field theories. 

We note that two reviews on quasinormal modes of 
astrophysical black holes and stars appeared in 1999 7,8 , 
and one recent review on quasinormal modes of black 
holes and branes^ appeared when we were writing this 
manuscript. When choosing the material for this re- 
view, we emphasized those questions which were not re- 
lated in detail in these reviews. Thus, we made a review 
of methods for separation of variables in the perturba- 
tion equations (Sec.|ll|, related presently used numerical 
and analytical methods for finding quasinormal modes 
(QNMs) fSec. lIII[) . and discussed in detail the stability of 
black holes (Sec. IVIII|) and late-time tails (Sec. IVII|) . In 
Sec. IIX| in addition to what has become the "standard" 
material on the hydrodynamic regime of field theories at 
strong coupling, we included a subsection that explains 
how ordinary WKB methods, usually used for QNMs of 
black holes, can be applied to finding the conductivity of 
holographic superconductors. 

In most cases, by "quasinormal modes" we mean "fre- 
quencies of quasinormal modes" and not the correspond- 
ing amplitudes. We will be mostly using the units 
h= c = G = 1. 



covariant form of the equation of motion is quite differ- 
ent for fields of different spin s in curved backgrounds. 
Thus, for a scalar field $ of mass /x in the background 
of the metric g^, the equation of motion is the general 
covariant Klein-Gordon equation 10 

(V"V V - M 2 )* = 0, 

where V„ is the covariant derivative. The above equation 
can be written explicitly as follows: 
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■■d v (g^V^gd^) - ^ = (s = 0). (2.1) 



For massive Dirac fields in a curved background g^ u , 
the equation of motion reads^i 

h a e a "(d,+T^+^ = 0, (a = ±1/2) (2.2) 

where \i is the mass of the Dirac field, and e a M is the 
tetrad field, defined by the metric g^ v : 



9nv — Vabe-ft e y , 



liv ab ll v 



a b 



cb a u 

°a: e e a 



(2.3) 



where r\ a \y is the Minkowskian metric, 7° are the Dirac 
matrices: 

{ 7 a , 7 b } = 2^ b , 
and is the spin connection^, 
1, 



r, 



(2.4) 



For massive vector perturbations we have the general 
covariant generalization of the Proca equations^. For a 
vector potential A M , one has 

VF^ - M % = 0, = - dvAp. (2.5) 

In a curved space-time these equations read 

1 



=d„( WV"V=ff) - ISA" = 0. (s = 1) (2, 



II. MASTER WAVE EQUATIONS 
A. Equations of motion 

From a theoretical point of view, perturbations of a 
black hole space-time can be performed in two ways: by 
adding fields to the black hole space-time or by perturb- 
ing the black hole metric (the background) itself. In the 
linear approximation, i. e., when a field does not backre- 
act on the background, the first type of perturbation is 
reduced to the propagation of fields in the background of 
a black hole, which is, in many cases, a general covari- 
ant equation of motion of the corresponding field. The 



When fj, = in the above form of the Proca equation, we 
obtain the Maxwell equation 



d v {(d a A a - d a A a )g a »g°»^-g) = 0. 



(2.7) 



There may be various generalizations of the massive 
scalar, spinor, and vector fields considered above. Thus, 
if we study perturbations of massive charged particles in 
scalar electrodynamics in a curved charged background, 
we have to deal with a complex scalar field 

(D V D V - /i 2 )* = 0, 

where D v — V v — ieA v is an "extended" covariant deriva- 
tive, and e is the charge of the particle. Finally, we find 
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that the equation of motion of the charged scalar held in 
a curved space-time readsi^ 

^(0 M * - ieA„*)) - ieA'Bv* 

-{[i 2 +e 2 A v A v )^ = 0. (2.8) 

In a similar fashion, the massive charged Dirac particle 
is described by the equation of motion with an extended 



derivative — » 



il a <{du +T U - ieA v ) + = 0. 



(2.9) 



Another type of perturbation, metric perturbations, 
can be written in the linear approximation in the form 



S T n 



9iw = 9% + Sg^v, 
T 9u 



D- 2 



2A 
D^2 l 



(2.10) 
?m- (2-11) 



Linear approximation means that in Eq. (|2.11[) the terms 
of order ~ Sg 2 and higher are neglected. The unper- 
turbed space-time given by the metric <7° is called the 
background. 



B. Separation of variables and radial equations of wavelike 
form 

The first step toward the analysis of the black hole 
perturbation equations is their reduction to the two- 
dimensional wavelike form with decoupled angular vari- 
ables. Once the variables are decoupled, a wave like 
equation for radial and time variables usually has the 
Schrodingcr-like form for stationary backgrounds, 



d 2 R 



+ V(r,w)R = u; 2 R 1 



(2.12) 



and can be treated by a number of sophisticated and 
well-developed numerical, analytical, and semianalytical 
methods. 

As a simple exercise, one can see that for the massless 
scalar field on the Schwarzschild background (<7° : gu = 

-g- 1 = 1 - 2M/r, g ee = gu sin" 2 6 = r 2 ), Eq. flU), 
after using a new variable 



dr 



1 - 2M/r 



The coordinate maps the semi-infinite region from 
the horizon to infinity into the (— oo,+oo) region and 
is, therefore, called the tortoise coordinate. 
The wave function 

*(t,r,6,ct>) = e- i " t Y e (6, C t>)R(r)/r, 

produces (|2.12l) with a potential 

^>4-^)(^ + ^f^). (-3) 



where s = 0, and I is the multipole quantum number, 
which arises from the separation of angular variables by 
expansion into spherical harmonics 

A e M6,(l>) = -e(t + l)Y i {6,<t>), 

exactly in the same way as it happens for the hydrogen 
atom problem in quantum mechanics when dealing with 
the Schrodinger equation. 

When s = 1 the effective potential (12.13[) corresponds 
to the Maxwell field. When s = 2 we obtain the effective 
potential of the gravitational perturbations of the axial 
type, which was derived by Regge and Wheeler ir>i^. 

The separation of variables, however, is not always 
so easy. The variables in perturbation equations can- 
not be decoupled for perturbations of an arbitrary met- 
ric. For this to happen, the metric must possess suf- 
ficient symmetry, expressed in existence of the Killing 
vectors, Killing tensors, and Killing- Yano tensor s 15 ' 16 . 
The choice of appropriate coordinates is crucial for sep- 
aration of variables^!. Having left readers with litera- 
ture on this subject^, in the next two subsections we 
considered some more generic and practically useful ex- 
amples of the separation of variables for fields of vari- 
ous spin s = 0, 1/2, 1,3/2, 2 in various black hole back- 
grounds. The two most efficient and general approaches 
to black hole perturbations will be related here: the gen- 
eral formalism of the spin s perturbations with the help of 
the Newman-Penrose tetrads^ and the gauge-invariant 
method for gravitational perturbations 2 ^. 



C. Separation of variables procedure for the 
Kerr-Newman-de Sitter geometry 

We considered the perturbation equations for fields of 
various spin in the the Kerr-Newman-(A)dS background, 
describing the charged, rotating black hole in asymptot- 
ically de Sitter (or anti-de Sitter) universe. This solu- 
tion described quite a wide class of black holes which 
possesses a number of parameters: the black hole mass, 
its charge, momentum, and the cosmological constant. 
Since the electromagnetic and gravitational fields couple 
through the electric charge of the black hole, there is no 
pure electromagnetic or pure gravitational perturbations 
unless the charge Q = 0, but one type of perturbation in- 
duces the other. With the help of the Newman-Penrose 
formalism shown below, Teukolsky in^i separated vari- 
ables in the perturbation equations for the Kerr geometry 
for massless fields of various spin. 

In Boyer-Lindquist coordinates the Kerr-Newman- 
(A)dS metric has the form, 



ds 2 



dr 2 

a: 



de 2 \ A e sin 2 e 

A7j~(l + a)V 



[adt-(r 2 + a 2 )dipf 



(1 + a) 2 p- 



■(dt — asm 2 9d(p) 2 , 



(2.14) 
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where 



A r = (r 2 + a 2 ) (l - ^r 2 ) - 2Mr + Q 2 (2.15) 

Oi 

= j(r - r+)(r - r_)(r -r + )(r -r_), 



A e = 1 + a cos 2 6>, 



p = r + tacost 



Aa 2 
3 ' 



Here A is the cosmological constant, M is the mass of 
the black hole, Q is its charge, and a is the rotation 
parameter. The electromagnetic field due to the charge 
of the black hole is given by 



A u daf = - 



Q> 



(1 + a) 2 p- 



■(dt ~ asm 2 9dtp). (2.16) 



In particular, we adopt the following vectors as the null 
tetrad: 



^ - 



(l + q)(r 2 +a 2 ) a(l + a) 

A ' ' 1 A 



^ ((1 + a)(r 2 + a 2 ), - A r , 0, a(l + a)) 



P\/2A g 



ia(l + a) sin 0, 0, Ag 



i(l + a) 



suit 



(2.17a) 



Because of the existence of the Killing vectors dt and d<j, 
for the unperturbed background, the time and azimuthal 
dependence of the fields has the form e -«M- m v) j the 
tetrad components of the derivative and the electromag- 
netic field are 



i"0u = 2? , 



„2 ^0! 



~^ fl £j, m^d^ = ^ £ , 



V2, 



V2p* 



(2.17b) 



/'M — — — 
M_ A r ' 
mM p = mM„ = 0, 



2p 2 



where 



i(l + a)i?f <9 r A r 

O,- 1 h n— : 



Making use of the above equations (|2.17l) and after 
some algebra the covariant equations of motion can be 
cast into the form with separated variables 



y/A e C\_ s VA e C s - 2(1 + a) (2s - l)awcos< 
-2a(s - l)(2s - 1) cos 2 9 + A] S s (6) = 0, 

A r DiX>] + 2(1 + a) (2s - l)iur - ^(s - l)(2s - 1) 



(2.18) 



-2(1 + a)eQKr + iseQrd r A r + e 2 Q 2 r 2 
A r 

-2iseQ - A] J2 s (r) =0, 



(2.19) 



where s is the spin of the field, and e is its charge. In the 
general case of the Kerr-Newman-de Sitter (and anti-de 
Sitter) space-times, spin and 1/2 fields can be separated 
for the above equations^ 2 -. For Kerr-de Sitter black holes 
(Q = 0), the separation is possible for fields of spin 0, 
1/2, 1, 3/2, and 2. 



D. Gravitational perturbations of D-dimensional black holes 

The Newman-Penrose tetrad approach to the separa- 
tion of variables given in the previous subsection was effi- 
cient for D = 4 black holes. The higher-dimensional cases 
were not decoupled until the recent papers of Ishibashi 
and Kodama^r— , although the D-dimensional general- 
izations of the Schwarzschild and Kerr solutions, called 
the Tangherlini 2 ^ and Myers-Perry^ solutions, respec- 
tively, have been known for a long time. 

For the purposes of string theory and higher- 
dimensional gravity mentioned in the introduction, we 
considered a general class of the higher-dimensional black 
holes, whose space-time can be represented as the direct 
sum of two-dimensional space-time (r, t) and (D — 2)- 
dimensional space M D ~ J\f 2 x KP~ 2 . The correspond- 
ing metric has the form 

ds 2 = g ab (y)dy a dy h + r 2 (y)daf D _ 2) , (2.20) 
where dajj_ 2 = r jij(z)dz' 1 "dz^ is a metric of the (D — 2)- 



dimensional complete Einstein space K. 



D-2 



Ri 



(D-3)K ll3 (K = 0,±1), 



and g a b is a static metric of the two-dimensional space- 
time TV 2 : 



C n = ^ + (1+ . Q)g +n^^ Sin 



a = d s - 



A e 
(l+a)H 



Ag sin 9 
fA~e sin( 



y/A~g sin 9 



with K = lo(t 2 + a 2 ) — am and H = —au> sin 9 - 



(2.17c) 



rn 
sin 9 



ab dy a dy b = -f(r)dt z + 



dr 2 
f(rY 



(2.21) 



In this subsection we assumed that the metric (|2.20l) is a 
solution of the vacuum Einstein equations with the cos- 
mological constant A. Therefore f(r) can be written in 
the form 2 ^ 



f(r) = K 



2M 



r fl - 3 (D-2)(I>-1) 



(2.22) 
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For K = 1, we have a regular black hole if 
2A 



(D-2)(D- 



_ M 2/(D-3) < (D_ 3 )/(iJ_l)(D-l 



-l)/(£>-3) 



For if = or —1, KP~ 2 may not be compact, and the 
space-time contains a regular black hole only for A < 0. 

We now start the analysis of linear perturbations of 
the above static D-dimensional space-times (12.20[) . The 
convenient perturbation variables are 

ipnv = - hg^/2, h^v = Sg^. 

Then, the perturbed Einstein equations can be written 
in the form 

-V 2 ^^2R liay p^+2V {p V a ^ u)a -V a V^ a pg tlu = 0. 

(2.23) 

These equations are invariant under the gauge transfor- 
mations 



-Vu6 - V„£„ 



(2.24) 



In order to extract the dynamics of the physical degrees 
of freedom from the above perturbation equations (|2.23[) 
the gauge freedom must be eliminated. 

The main technical challenge of the separation of vari- 
ables problem for gravitational perturbations is that 
Eq. (|2.23p is a set of coupled equations with D(D + l)/2 
unknown functions. For any metric that can be written 
in the form (|2.20[) . the separation of variables is possible 
with the help of the following tensorial decomposition of 
h^ v . According to the tensorial behavior on KP~ 2 , h^ u 
can be divided into the following components h a b, h a i 
and hij. Here letters a, b, and c are used for the radial 
and time coordinates, while coordinates on KP~ 2 are de- 
noted by i, j .. . Then, the vector and tensor components 
can be decomposed as follows^ 



iTij 



(2.25) 
(2.26) 



h al = A*a + *ai! D l ^! ai = 0, 



/) >!'/ , = 0, l)-A> i , 



0. 



(2.27) 



where Di is the covariant derivative with respect to the 
metric jij on KP~ 2 , A = D ■ D. In this way one can 
obtain the following three groups of variables, which are 
classified by their transformation law with respect to the 
coordinate transformations on KP~ 2 : 



• scalar-type variables ^ a b 

• vector-type variables 



* Q , and *t; 
and * Ti ; 
• and a tensor- type variable ^nj- 



The Einstein equations written in this form split into 
three subsets each of which contains only variables be- 
longing to one of the above three sets of variables. Then, 
the gauge-invariant variables can be constructed from 



the harmonic expansion coefficients. Through this pro- 
cedure the Einstein equations can be transformed to a 
set of gauge-invariant equations with a small number of 
unknown functions^. In particular, for tensor perturba- 
tions, this procedure gives a single wave equation for a 
single unknown function (|2.12[) . 

When D — 4 we have a particular case of the gen- 
eral approach. The scalar type of gravitational perturba- 
tion was called "polar" by Chandrasekhar™. The polar 
(scalar) perturbations are responsible for deformations of 
the black hole horizon. The vector type is known as "ax- 
ial" in four dimensions. This type of perturbation is con- 
nected with small rotations of the black hole. The tensor 
perturbations are not dynamical and can be gauged off in 
four dimensions. Note that it is important not to confuse 
the scalar type of gravitational perturbations (sometimes 
simply called "scalar perturbations" ) with perturbations 
of the test scalar field (|2.ip . 



III. METHODS FOR QUASINORMAL MODES 
CALCULATIONS 

A. Definitions and main properties of the quasinormal 
modes 

Quasinormal modes are solutions of the wave equa- 
tion (|2.12p , satisfying specific boundary conditions at the 
black hole horizon and far from the black hole. At the 
event horizon this boundary condition is a requirement 
of the pure ingoing waves 



pure ingoing wave, 



(3.1) 



Another boundary condition, imposed at spatial infinity 
or at the de Sitter horizon for asymptotically de Sitter 
black holes, is different in the astrophysical and string 
theory contexts. For astrophysical purposes one has 



^ ~ pure outgoing wave, r* — > +oo. 



(3.2) 



In string theory the boundary condition at infinity de- 
pends on the perturbed field under consideration. For 
the simplest scalar field the wave function ^> must van- 
ish at infinity, while for higher spin fields some gauge- 
invariant combination of field components, dictated by 
the AdS/CFT correspondence^, must vanish. Thus in 
string theory the boundary condition at infinity is the 
Dirichlet one 

*^0, r^oo (s = 0), (3.3) 
gauge inv. comb. — > 0, r — > oo (higher s). 

Remembering that ^ ~ e" ,ul , we write the QNM fre- 
quencies in the following form: 



UJ = LORe - IWlr, 



(3.4) 



Here u>R e is the real oscillation frequency of the mode 
and uJim is proportional to its damping rate. Positive 
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LJi m means that ^ is damped, negative w/ m means an 
instability. 

We now review the main properties of the quasinormal 
modes of black holes. 

1. Quasinormal modes of non-AdS black holes do not 
form a complete set. Therefore, the signal can- 
not be represented as a sum over the quasinormal 
modes all of the time. Thus, at asymptotically late 
time, t — ¥ oo, the power law or exponential tails 
usually dominate in a signal (see Fig. [3]). An ex- 
ception to this rule is the asymptotically AdS black 
holes, for which, due to their special boundary con- 
ditions, quasinormal modes dominate at all time 3 ^. 

2. The quasinormal frequencies do not depend on a 
way by which the black hole or a field around it 
was perturbed. Thus, quasinormal modes are com- 
pletely determined by a black hole's parameters 
and were called, therefore, "fingerprints" or "foot- 
prints" of black holes. 

3. For Kerr black holes as well as for other astrophys- 
ical or string theory motivated cases, quasinormal 
modes form a countable set of discrete frequencies 

4. Quasinormal modes calculated in the linear approx- 
imation are in good agreement with those obtained 
by the fully nonlinear integration of the Einstein 
equations, at least at sufficiently late tim e 31 ' 32 . 

Practically, it is important to calculate QNMs with a 
high accuracy because considerable changing of black hole 
parameters frequently changes quasinormal frequency 
just by a few percent. Therefore, numerical and semian- 
alytical methods for solving the quasinormal eigenvalue 
problem with high accuracy have gained considerable at- 
tention during the past years. 



such that the potential is invariant under these transfor- 
mations 



V(r*,a) = V(—ir*,a'). 



(3.7) 



The wave function ^ and the quasinormal frequency u> 
transform as 

tf(r„, a) = $(-tr»,a'), w(a) = Sl(a') (3.8) 

Then, the wave equation for $ and the boundary con- 
ditions will read 

^+(-n 2 +V)<f> = 0, $ ~ e^ nr * , r* too. (3.9) 

These boundary conditions correspond to a vanishing 
wave function at the boundaries, so that the QNM prob- 
lem is now reduced to the bound states problem for an in- 
verse potential V — > —V, which is smooth potential gap, 
approaching some constant values at the infinite bound- 
aries. This gap can be approximated by a potential for 
which we know the analytic solution of the wave equa- 
tion: It is the Poschl-Teller potential 33 , 



Vpt 



V 



cosh a(r* — r°) 



(3.10) 



Here Vb is the height of the effective potential and 
—2Voa 2 is the curvature of the potential at its maximum. 
The bound states of the Poschl- Teller potential are well 
known 34 



1 



1 , v 



4 (a') 2 



1/2N 



(3.11) 



The quasinormal modes uj can be obtained from the in- 
verse transformation a' = ia, 



B. The Mashhoon method: Approximation with the 
Poschl-Teller potential 

Here we relate possibly the easiest method for calcula- 
tion of the quasinormal modes, which is also illustrative 
as to the physical essence of the problem. This method 
was suggested by B. Mashhoon^ 3 .. 

We start from the usual wavelike equation, with an 
effective potential, which depends on some parameter a: 



— + {u 2 -V{r^am 



0. 



(3.5) 



According to (|3.1[ I3.2j) the astrophysically relevant quasi- 
normal boundary conditions for asymptotically flat black 
holes are 



— > zboO. 



(3.6) 



Because of "symmetric" boundary conditions for the 
QNM problem at both infinities r» — > ±oo, it is reason- 
able to consider transformations r* — > — ir* and p — > p' , 



±\V 



1 



ia n + 



1 



0,1,2,. 



(3.12) 

Technically one has to fit a given black hole potential 
to the inverted Poschl-Teller potential. In the forthcom- 
ing sections we will see that for the low-lying QNMs, in 
the majority of cases the behavior of the effective poten- 
tial is essential only in some region near the black hole, so 
that the fit of the height of the effective potential and of 
its curvature is indeed sufficient for estimation of quasi- 
normal frequencies. 

This method gives quite accurate results for the regime 
of high multipole numbers t, i.e. for the eikonal (geo- 
metrical optics) approximation. In particular, for gravi- 
tational perturbations of the D — 4 Schwarzschild black 
holes, the fundamental quasinormal modes obtained by 
the Mashhoon method have a relative error of not more 
than 2% for the lowest multipole £ — 2, and of about 
fractions of a percent for higher multipoles. 

There are cases when the effective potential of a black 
hole is exactly the Poschl-Teller potential. These are 
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Schwarzschild-de Sitter and Reissner-Nordstrom-de Sit- 
ter black holes with extremal value of the A-ter m 35 ' 36 . 
After some modifications the Mashhoon method can be 
used for rotating black holes^I. Other exactly solvable 
potentials and their applications to the analysis of the 
QN spectra were reviewed in 3 ^. 



C. Chandrasekhar-Detweiler and shooting methods 

By the substitution 

.-«,(< /%*.), (3.13) 

the wave equation (|2.12l) can be reduced to the Riccati 
equation 

id<f>/dr* + u? - $ 2 - 7(r*) = 0. (3.14) 

The quasinormal boundary conditions (|3.6|) give 

$ -> +cj, r* -> +oo, (3.15) 
$ -> -w, -> -oo. (3.16) 

Then in order to obtain QNMs one needs to integrate 
the Riccati equation numerically 3 ^. The easiest way is to 
"shoot" the two asymptotic solutions at the horizon and 
in the far region at some common intermediate point, 
which is usually the peak of the potential barrier. The 
shooting procedure is the following. We take some fixed 
(complex) ui with if? oc e 1 " 1 ** and integrate from infin- 
ity backward to some intermediate value of r* = r*. 
Then, for the same u we integrate the wave equation 
with if? oc e~~ lur * from the horizon until the chosen inter- 
mediate value of rjjj. In the common point, the Wronskian 
of the two solutions must vanish, which gives the equa- 
tion for the quasinormal modes 3 ^. This shooting scheme, 
when applied to the Riccati equation, works quite well 
for lower modes, yet, there is a technical difficulty in this 
procedure when one uses it for the usual second-order 
wave equation (|2.12|) . For large positive r*, the solu- 
tion if? with asymptotic $ oc e lur * gets contaminated by 
the admixture of waves oc e~ lulr * and vice versa; for large 
negative r* the solution with asymptotic if? oc e~ lur * gets 
admixture with waves oc e^ zwr * . 

Since for the AdS case the asymptotical behavior of 
the function is not exponential it is possible to adopt the 
shooting method for the second-order differential equa- 
tion. Namely, we integrate Eq. (|2.12[> . starting from the 
event horizon, where the quasinormal boundary condi- 
tion is applied for a fixed value of uj. At large value of r 
we find a fit of the function if? by the two linearly inde- 
pendent solutions, in which the series expansions can be 
found analytically up to any order at spatial infinity. One 
of these solutions if! d vanishes at infinity, satisfying the 
boundary condition inspired by string theory, while the 
other solution 'Fjv grows. The fitting procedure allows us 




x 



FIG. 1 The three regions separated by the two turning points 
Q{x) = 0. 

to find the contributions for each of these solutions for 
the particular value of uj: 

*(r) « C D (uj)if? D {r) + C N (uj)^ N (r), r > r+. 

The requirement of vanishing of the function if? at spa- 
tial infinity implies that C/v(w) = 0, if the frequency ui is 
a quasinormal mode. Because of growing of the function 
if 1 n{t), the fitting procedure provides good accuracy for 
the corresponding coefficient and allows us to find accu- 
rate values of quasinormal modes by minimizing Cn(uj). 

For the asymptotically flat and de Sitter cases the 
exponentially growing solution satisfies the quasinormal 
boundary conditions, while the other solution decays ex- 
ponentially for large r. The exponentially decaying con- 
tribution is difficult to calculate. That is why the shoot- 
ing method for the second-order equation does not allow 
us to find quasinormal modes for stable non-AdS back- 
grounds. 

Yet, if the quasinormal spectrum contains a growing 
mode (uji m < 0), the corresponding eigenfunction if? 
must be zero at the spatial infinity. Therefore, we are 
able to use the shooting method to check stability of 
black holes. In order to do this one can check for all 
possible values of w, for which the corresponding if? does 
not vanish at spatial infinity (or at the cosmological hori- 
zon). 

In practice, we consider only the parametric region of 
uj, where we were unable to prove analytically that the 
growing mode is not consistent with the equation and 
the quasinormal boundary conditions. In this region we 
calculate the prefactor for the growing solution. If the 
pref actor is zero for some value of uj, this is a growing 
quasinormal mode. 

D. WKB method 

Here we describe the semianalytic technique for finding 
low-lying quasinormal modes, based on the JWKB ap- 
proximation. The method was first applied by Schutz and 
Will to the problem of scattering around black holes^. 
This method can be used for effective potentials which 
have the form of a potential barrier and take constant val- 
ues at the event horizon and spatial infinity (see Fig. [1} . 
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The method is based on matching of the asymptotic 
WKB solutions at spatial infinity and the event horizon 
with the Taylor expansion near the top of the potential 
barrier through the two turning points. 

We rewrite the wave equation (|2.12[) in a slightly dif- 
ferent form 



dx 



— + Q(x)*(x)=Q. 



(3.17) 



Here we used the new designations x = r* and Q{x) = 
uj 2 — V. The asymptotic WKB expansion at both infini- 
ties has the following general form: 



exp 



(3.18) 



We introduced the WKB parameter e in order to track 
orders of the WKB expansion. Substituting the expan- 
sion (|3.18p into the wave equation p. 171) and equating 
the same powers of e, we find 

S (x) =±i J Q(v) 1/2 dv, (3.19) 



Si(a;) = --hiQ(z) 



(3.20) 



The two choices of sign in (|3 . 19[) correspond to either 
incoming or outgoing waves at either of the infinities 
x = ±oo. Thus, when x — > +oo (region I of Fig. [TJ, 
Q(x) — > uj 2 in the dominant order, so that So — > +icux 
for the outgoing to the infinity wave and So —itox for 
the incoming from infinity wave. In a similar fashion, at 
the event horizon x — > — oo (region III), So — >• +iuix is for 
wave incoming from the event horizon, while So — > —itux 
is for a wave outgoing to the event horizon. We designate 
these four solutions as 'JK , ^SfL, and ^ I J I respec- 

tively for plus and minus signs in So in regions I and III. 
Thus, 



-oo, 



(3.21) 



If the turning points are closely spaced, i.e. if 
—Q(x)max <C Q(±oo), then the solution in region II can 
be well approximated by the Taylor series 

Q(x) = Qo + \Q'a(x - x ) 2 + 0((x - x ) 3 ), (3.26) 

where xo is the point of maximum of the function Q{x), 
Qo — Q{ x o), an d Qo is the second derivative with respect 
to x at the point x = xq- Region II corresponds to 



I I , /— 2Qo _ 1/2 



region II. (3.27) 



The latter relation gives also the region of validity of the 
WKB approximation: e must be small. 
We introduce new functions 



k= l -Ql t={Ak) l ' i e™/\x-x Q ), 



1 



/Ml/2 



Then, the wave equation (|3 . 1 T[) takes the form 
d 2 * / 11 



dt 2 



2 4 



t z * = 0. 



(3.28) 
(3.29) 

(3.30) 



The general solution of this equation can be expressed in 
terms of parabolic cylinder functions D v {t), 

* = AD v (t) + BD- y -i(it). (3.31) 

Large \t\ asymptotics of this solution are 

-3i,ir(„+l) </+! S-IV+lI 'fc 1/2 ("-^0) 2 

+(A + B(27r) 1 / 2 « _ "" r/a /r(i/ + l))x 



e^(4/c)T(a; - x ) l/ e ""^ X0> , a;>X2, (3.32) 



,J7J 



,iu 



The general solutions in regions I and III are 

region I, 



(3.22) 9 ~ A e - 3 ^/ 4 (4fc)^ 4 (x - Xo )" e -^ 1/2 (*-°) 2 /2 + 

/B-iA(2 7 r) 1 / 2 e- WTr / 2 \ *.wi> . j ,, -wd 

(3.23) J e^(4fc)^" x 



'-'out — > 



region 



III. 



(3.24) 



The amplitudes at +oo are connected with the ampli- 
tudes at — oo through the linear matrix relation 



7 iii 

'-'out 



Sn S12 

S21 S22 



7 1 

''out 

zL 



(3.25) 



Now we need to match both WKB solutions of the form 
(|3. 18|) in regions I and III with a solution in region II 
through the two turning points Q(x) = 0. 



(x-x r { " +1) e* kl/2(x - xo)2/2 , x< Xl (3.33) 



Equating the corresponding coefficients in (|3.32l) . 
p.33|) and eliminating A and B, we obtain the elements 
of S matrix, 



7 iii 

'-'out 

rzlll 



R- 2 (2^) 1 ' 2 

r(-») 



t R 2 e""(2^) 1 ' 2 




(3.34) 
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where 



E= ( I/ + l)(-+l/2)/2 e -(v+l/ 



2)/2 



(3.35) 



When expanding to higher WKB orders, S matrix has 
the same general form (|3.34p . though with modified ex- 
pression for R, which still depends only on v. We note 
that for a black hole there are no waves "reflected by 
the horizon" , so that Z^ 1 = 0, and due to quasinormal 
mode boundary conditions, there are no waves coming 
from infinity, i. e., Z\ n = 0. Both these conditions are 
satisfied by (|3.34j) . only if 



r(— v) = oo, 



(3.36) 



and, consequently, v must be an integer. Then, from the 
relation (|3.29l) we find 



n+ l - = -iQ /{2Q'tf/ 2 , 



n = 0,1,2,.. 



(3.37) 



The latter relation gives us the complex quasinormal 
modes labeled by an overtone number n at the first WKB 
order—. Later this approach was extended to the third 
WKB order beyond the eikonal approximation by Iyer 
and Will^i and to the sixth order by Konoply a 42 ' 43 . In 
order to extend the WKB formula to higher orders, it is 
sufficient to take higher orders in the e WKB series (|3.18[) 
and to take an appropriate number of consequent terms 
in the Taylor expansion (|3.26j) . Since the S- matrix (|3.25j) 
depends only on its elements Sij can be found simply 
by solving the interior problem in region II at higher or- 
ders in eQ, and without explicit matching of the interior 
solution with WKB solutions in regions I and III at all 
the same orders. 

Going over from Q to the effective potential V, the 
sixth order WKB formula has the form 4 ^ 



^ 2 - Vb) 



6 

E 

i=2 



A s: -=n + i, n = 0,l,2,... (3.38) 



and the reflection coefficient is R = 1 — T. 

It was shown i n 42 i 43 that the WKB formula, extended 
to the sixth order, gives the relative error about 2 or- 
ders less than that of the third WKB order. Strictly 
speaking the WKB series converge only asymptotically, 
so that the consequent decreasing of the relative er- 
ror in each WKB order is not guaranteed. Therefore 
it is reasonable to develop a modified WKB technique 
in the so-called optimal order 45 . The latter gives bet- 
ter results for moderately higher overtones n and espe- 
cially when n > I. Yet, in many cases when n < £ 
the usual sixth order WKB formula gives much better 
results than the optimal order approach. For exam- 
ple, for the fundamental (odd parity) gravitational mode 
of the Schwarzschild black hole, the sixth WKB QNM 
ojM = 0.3736 — 0.0890i is much closer to the accurate 
numerical value ujM — 0.3737 — 0.0890i than the optimal 
(third) order value ujM = 0.3734 - 0.089H given by the 
phase-integral approach^ 5 .. This is in agreement with the 
general experience of various applications of the WKB 
method: Extension to higher orders frequently gives bet- 
ter accuracy than expected 4 ^. 

In some cases the WKB approach related here needs 
modifications. For instance, when considering a massive 
scalar field in a black hole background, the effective po- 
tential has a local minimum far from a black hole. This 
local minimum induces two changes in the WKB proce- 
dure: First, there are three turning points which sepa- 
rate all space into four regions, so that three matchings 
are required 4 ^. Second, an influent subdominant term in 
the asymptotic WKB expansion at spatial infinity (|3.24[) 

48 49 

appear s ' . 

The higher-order WKB approach proved to be use- 
ful for finding lower overtones of the quasinormal 
spectrum 5 ^— and is in good agreement with accurate 
numerical data as was shown i n 54 ' 55 . The higher-order 
WKB corrections lead also to better accuracy for the re- 
flection and transmission coefficients 56 . 



where the correction terms depend on the value of 
the effective potential and its derivatives (up to the i— th 
order) in the maximum. The explicit form of the WKB 
corrections can be found in^i (A2, A3) and in^i (A4, A5, 
A 6 ). 

In addition to solving the quasinormal mode problem, 
the S'-matrix allows us to solve the standard scattering 
problem, which describes tunneling of waves and particles 
through the potential barrier of a black hole. One can 
easily check from Eq. (|3.34[) that for real Q(x) (for real 
energy of the incident wave and spherically symmetric 
backgrounds) 



Sn — S22, S11 — S%i, IS21I 



The transmission coefficient is 



1. 



(3.39) 



E. Integration of the wavelike equations 

We rewrite the wavelike equation (|2. 12|) without im- 
plying the stationary ansatz ~ e~ lult ) 



<9 2 $ 



d 2 <S> 

dx 2 



(3.41) 



where x is the tortoise coordinate. The technique of in- 
tegration of the above wave equation in the time do- 
main was developed by Gundlach, Price and Pullin 57 . 
We shall rewrite the wavelike equation p. 411) in terms 
of the so called light-cone coordinates du — dt — dx and 
dv = dt + dx: 



7lII\2 



T = 



q-l 
D 11 ! 



(3.40) 



dudv 



+ V(u,v) $(u,v) = 0. 



(3.42) 
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FIG. 2 The integration grid. Each cell of the grid represents 
an integration step. The thick points illustrate the choice of 
(S, W , E, and N) for the particular step of the integration. 
The initial data are specified on the left and bottom sides of 
the rhombus. 



In these coordinates the operator of the time evolution is 



exp h 



d_ 

du 



exp h 



T ( exp ( h lL 



d_ 

dv 

exp h 



- 1 + 
d \\ d 2 



dv J J dudv 



0(h 4 ). 



Acting by this operator on $ and taking account of 
(|3~I2j) . one finds 

$(jV) = $(W) + $(E) - $(S) - 

^S)($(WO + $(i;))+0(/i 4 ), (3.43) 

where we introduced letters to mark the points as follows: 
S = (u, v), W = (it + h, v ), E = (u, v + h), and N = 
(u + /i, u + h). Equation Q3.43P allows us to calculate 
the values of $ inside the rhombus, which is built on 
the two null-surfaces u = uq and v — vq (see Fig. [5]), 
starting from the initial data specified on them. As a 
result we can find the time profile data {$(£ = to), 3>(f = 
<o + h), <J>(t = to + 2h), . . .} in each point of the rhombus. 
These values of the function can be used for calculations 
of quasinormal modes. 

In this way, one can obtain a time-domain profile of 
the perturbation, an example of which is shown in Fig. [3] 
When looking at Fig. [31 a natural question is to which 
value of the quasinormal frequency the above profile cor- 
responds. The answer can be found with the help of the 
Prony method of fitting of the profile data by superposi- 
tion of damped exponents (see, e. g-r^) 



$(t)~53Ci, 



(3.44) 




250 



FIG. 3 An example of a time-domain profile for the 
Schwarzschild black hole gravitational perturbations (1 = 2 
vector type, in the point r = llr+). 



We suppose that the quasinormal ringing epoch starts at 
to = and ends at t = Nh, where TV is an integer and 
N > 2p — 1. Then the formula (|3.44j) is valid for each 
value from the profile data: 

p p 
x n = <S>(nh) = Cje' l ^ nh = C 3 z™. (3.45) 

The Prony method allowed us to find zi in terms of known 
x n and, since h was also known, to calculate the quasi- 
normal frequencies Wj. In order to do this, we defined a 
polynomial function A(z) as 



A(z) = l[(z-z j ) = 



p 

E< 

m=0 



a = 1. (3.46) 



Consider the following sum: 

p p 



p 



m^n—m — / , / J j 

m— m— J=l 



P 



j — l m— 

Since ao = 1 we find 

p 



E 

rn — 1 



m-^n — m 



(3.47) 



Substituting n = p, N into Eq. (|3.47p we obtain N 
p + 1 > p linear equations for p unknown coefficients a n 
We rewrite these equations in the matrix form 



Xp_i Xp—2 
Xp Xp — x 



2^0 \ 

Xi 



a 2 



( x P \ 

Xp+l 



\ x N ^i x N ^ 2 ■ ■ ■ XN-p J \a p J \ x N J 
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Such a matrix equation 

Xa = —x 

can be solved in the least-squares sense 
a = -(X+X)- 1 X+x, 



V(r) 



(3.48) 



where X + denotes the Hermitian transposition of the 
matrix X. 

Since the coefficients a m of the polynomial function 
A(z) are found, we can calculate numerically the roots 
Zj of the polynomial and the quasinormal frequencies u)j : 



UJj 



While the above described integration scheme (|3.43|) 
was efficient for asymptotically flat or de-Sitter black 
holes, for asymptotically AdS black holes its convergence 
is too slow. For the latter case, an alternative integration 
scheme^, 
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[V(S)*(S) + V(E)¥(E) + V(W)V(W)} , (3.49) 



is more stable. 



F. Fit and interpolation approaches 

For asymptotically flat (or de Sitter) black holes, the 
wavelike equation usually has the form (12.121) . The gen- 
eral solution of Eq. (|2. 12|) at infinity can be written in 
the form 



* = A in tp in + A out ijj 



(3.50) 



The quasinormal modes, by definition, are the poles of 
the reflection coefficient A out /Ai n . In astrophysical con- 
text one does not work with the true asymptotic spatial 
infinity r* = oo, but a more practical notion of the "far 
zone" can be used instead. The far zone, or "astrophysi- 
cal infinity", is separated from a black hole by a distance 
that is much larger than the black hole radius r ^> r + . 
The basic intuitive physical idea of what we call here 
the "fit and interpolation approaches" is that scattering 
properties as well as quasinormal modes of black holes (in 
an astrophysical context) depend only on the behavior of 
the effective potential in some region near the black hole 
and do not depend on behavior of the potential at infin- 
ity— The immediate confirmation of this idea would 
is the WKB formula (|3. 381) or the Mashhoon formula for 
w that depends only on derivatives of the effective po- 
tential in its maximum, i. e., where the main process of 
scattering is localized. However, such a confirmation is 
almost trivial as it is based on the Taylor expansion of 
the effective potential near its maximum, which takes ac- 
count of the behavior of the potential only near its peak. 




FIG. 4 Potential for electromagnetic perturbations near the 
Schwarzschild black hole (r+ = 1, £ = 2) and the same poten- 
tial interpolated numerically near its maximum. Despite the 
behavior of the two potentials being different in the full region 
of r, except for a small region near the black hole, low- lying 
quasinormal modes for both potentials are close. 



Nevertheless, one can check that the above statement 
about the dominance of the "near zone" in scattering 
is true by considering the well-known potential for the 
Schwarzschild black hole V(r) and also two other poten- 
tials which lay perfectly close to the Schwarzschild po- 
tential near its maximum, but have strikingly different 
behavior far from the black hole. Then it can be shown 
that all three potentials, quite surprisingly, produce the 
same quasinormal modes. 

The two alternatives to the Schwarzschild potential arc 
chosen in the following way. We plotted the function 
V(r) of the well-known analytical Schwarzschild potential 
(I2.13[) and found numerical values of plot points of V(r) 
near its peak, which served as a numerical data for our 
first alternative potential. We call it Vint-, which is an 
interpolation of these points near the potential maximum 
by cubic splines (see Fig.[3J. The second potential, called 
Vfu, is a fit of the above plot's points near the maximum 
by a ratio of polynomial functions. From Fig. [4] one can 
see that far from the black hole, Vint is indeed different 
from the Schwarzschild potential because it diverges at 
large r. 

The method that is not "trivial" for our purpose is 
the time-domain integration described in Sec. IIII.E1 In- 
tegrating the wave equation with the above two effective 
potentials in time domain and comparing the obtained 
QNMs with their WKB values one can show that in- 
deed the quasinormal modes of all three potentials are 
the same up to the small numerical error—. 

The advantage of this approach is that one does not 
need to know the analytic form of the effective potential 
(or even of the background metric). Thus, it can be ap- 
plied to a number of problems of astrophysical interest, 
where a black hole is surrounded by some distributions 
of matter or, for instance, for some alternative theories of 
gravity^ 3 , where an exact solution for a black hole cannot 
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be found. What one needs is only a numerical solution 
for the background metric near a black hole. Strictly 
speaking, the notion "near" depends on the distance at 
which the black hole effective potential reaches its maxi- 
mum. Usually it is sufficient to know a numerical solution 
for the black hole metric in the region that starts from 
the event horizon and finishes at about 4—5 black hole 
radiiffi. 

G. Frobenius method 

Eq. (|2.12[) can be written in a slightly different form 

(^2+P^Tr +q{r) ) R{r)= ^ (3 - 51) 

where the functions p(r) and q(r) depend on the eigen- 
frequency uj. 



We started from the analysis of singularities of 
Eq. (|3.51[) . There are two points which are always sin- 
gular: the event horizon r = ?*+ and the cosmological 
horizon (or spatial infinity) r = r^. Usually, there are 
also other singular points, which depend on p(r) and 
q{r). By definition, quasinormal modes are eigenvalues 
u) which satisfy the boundary conditions corresponding 
to the outgoing wave at spatial infinity and the ingoing 
wave at the horizon. Thus, we are able to determine the 
function ^(r) as a multiplication of a divergent (at these 
points) function by a series, which is convergent in the 
region r + < r < r^. If p(r) and q(r) are rational func- 
tions of r, we can construct such a series in terms of the 
rational functions: 



R(r) 



r 



= < 



r-r 



e lUr (r~r Q ) c 



ro 

r+ 
■ r 



—ia oo 

E 

—ia oo 



E & * 

k=0 



ro r 



r+ 
■ r 



< oo, 



oo. 



(3.52) 



f2, a, and a are defined in order to satisfy Eq. ()3.51j) in 
the singular points r = r + and r = . The quasinormal 
boundary conditions fix Re(f2) and Re(a), which must be 
chosen to be the same sign as Re(w). 
The Frobenius series is 



u(z) 



OO 

E 

fc=0 



b k z k , 



r - r + - r 



r -r r c 



r+ 



(3.53) 



If all the singular points of Eq. (|3.51|) satisfy \z\ > 1, the 
series (|3.53[) is convergent at z = 1 (r = r^) if and only if 
the value of ui is the eigenfrequency of Eq. (|3.51|) . If there 
is at least one singular point inside the unit circle, one 
has to continue the Frobenius series (|3.53[) through some 
midpoints (see Sec. lIII.j]) . in order to test the convergence 
at spatial infinity or the cosmological horizon. 

Note that the definition of z contains an arbitrary pa- 
rameter r < r + . In most cases, it can be chosen in order 
to move all of the singularities outside the unit circle. 



H. Method of continued fractions 



Substituting (|3.52[) into (|3.51[) . one can obtain an TV- 
term recurrence relation for the coefficients bi. 

min(N—l,i) 

J2 <$P(u)&i-i=0, fori>0, (3.54) 

3=0 



where the coefficients Cjp(w) (0 < j < min(N — l,i)) 
depend on lo. 

We now decrease the number of terms in the recurrence 
relation 



min(k.i) 

£ cg +1 V)^-=0 

3=0 



(3.55) 



by one, i. e., we find (uj) which satisfy the following 
equation 

min(k-l.i) 

cJJ = 0. (3.56) 

3=0 

For i > k we can rewrite the above expression as 



( k ) , , 



E 4-i,-i (w)6i-i=0. (3.57) 



Subtracting (13.57)) from (|3.55l) we find the relation (|3.56[) 
explicitly. Thus, we obtain 

cfj{u) = c^ +1) (a;), for j = 0, or i < k, 

O+i)/ \ O) / A 
(fe)/ \ 0+1)/ n C k. t ( w ) c J-l,i-l( w ) 
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This procedure is called the Gaussian eliminations and 
allows us to determine the coefficients in the three-term 
recurrence relation 



c ol b i + c m h-i + 41 & t -2 =0, i > 1, (3.58a) 



,(3) 



Co I &i + ci 3 ! 6 = 0, 



(3.58b) 



numerically for a given u> up to any finite i. The com- 
plexity of the procedure is linear with respect to i and 
N. 

If the Frobenius series is convergent, we are able to find 
6i/6 from Eq. (|3.58b|) and substitute it into Eq. (|3.58aj) 

(3) (3) (3) (3) (3) (3) 



C 2,2 C 0,2 C 2,3 C 0,3 C 2,4 



61 J 

b J 3 ) r ( 3 )_ r (3)„ J3) 

C n 1 2 1.3 c l,4 



u 0,l 



Finally we find 

J3)J3) (3) (3) 
n _ (3) c 0,l c 2,2 C 0,2 L 2,3 

c l,2 c l,3 

The latter relation can be inverted n times to give 



(3.59) 



.(3) (3) 



(3) C 2,n C 0,ra-l C 2,n-l u 0,«-2 



.(3) „(3) 
1 



-1,71+1 



c (3) - 

c l,n-l 

J3) „0) 



„(3) _ 
-l,n-2 

.(3) J3) 



.(3) (3) 
"2.2 L 0.1 

1,1 



,(3) 
-l.n+2 



J3) 
-l,n+3 



(3.60) 



Equation (13.60)) with an infinite continued fraction on the 
right-hand side can be solved numerically by minimizing 
the absolute value of the difference between the left- and 
right-hand sides. The equation has an infinite number of 
roots (corresponding to the QN spectrum), but for each 
n, the most stable root is different. In general, we have 
to use the n times inverted equation to find the n-th 
QN mode. The requirement that the continued fraction 
is convergent itself allows us to truncate its length by 
some large value, always ensuring that an increase in this 
value does not change the final results within a desired 



precisio: 



,64 



I. IMollert improvement 

It turns out, that the convergence of the infinite con- 
tinued fraction becomes worse if the imaginary part of to 
increases with respect to the real part. This means that 
in order to calculate higher overtones, we must increase 
the depth of the continued fraction, which dramatically 
increases the time of computation. The convergence is 
poor also if ro in (I3.52|) is not a singular point. This fix- 
ing of ro is necessary to move all the singular points for 
higher-dimensional Schwarzschild black holes outside the 
unit circle \z\ < 1. 



The problem of slow convergence was circumvented 
ir>££ for the three-term recurrence relation and general- 
ized for higher N We consider the relation 



h c (3) c (3) c (3) 

On D C 2,n+1 c 0,n+l c 2,n+2 

= tin — 



6n-l 



,(3) 



,(3) 



^l.n+l u l,n+2 

which can be expanded for large n as 

i?„H^CoH + ^M + ^M 

y/n n 



(3.61) 



(3.62) 



In order to find the coefficients Cj of (|3.62p we di- 
vide Eq. p.54j) by 6,j_7v+i an d use the definition R n — 
— 6 n /6 n _i- Then one finds the equation with respect to 



N-l N-2-j 
0=0 fc=0 



(3.63) 



For large n we have c^ n (u>) oc n , thus, substituting the 
expansion (|3.62p into (|3.63p . we find 

N-l 

lim iEHNJMC 1 "''^) = 0- (3-64) 

j=0 



In the general case Eq. (|3.64[) has N — 1 roots (in fact, 
there are repeated roots among them) . If the series has 
a unit radius of convergence, one of the roots is always 
Cq = — 1 (it is also a repeated root) (|3.52[) . Other roots 
appear due to the existence of additional singular points 
of Eq. (|3.51j) . Thus, we choose Co = -1. 

After fixing Co = — 1 one can find an equation with 
respect to C\. In order to fix the sign of C\ we can use 
the convergence of the series (|3.52j) at z = 1. Therefore, 

lim b n = 0, i. e. $N: Vn > N, |6„| > |6„_i|. 
Since for large n 

& » t> r Ci _ Ci 

7 —tin ~ — <^0 7= — 1 1= , 

bn-x \Jn y/n 

we find out that the real part of C\ cannot be negative. 

Once the sign of C\ is fixed, the other coefficients in 
P-62|) can be found from Eq. f|3.63[) step by step without 
encountering indeterminations. 

As one calculated the coefficients Cj, the expansion 
(I3.62p could be used as an initial approximation for the 
"remaining" infinite continued fraction. In order to en- 
sure the convergence of p.62p for a given value of u, one 
has to start from the found approximation deeply enough 
inside the continued fraction (|3 . 60[) . The expansion gives 
a good approximation for R n . Therefore, the required 
depth is less than it would be if we had started from 
some arbitrary value. 
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J. Continuation of the Frobenius series through midpoints 

Consider now the case when one cannot fix the param- 
eter r in (|3.52l) in such a way that all the singularities, 
except r = r + and r — r^, can be removed outside the 
unit circle \z\ < 1. In other words, there is at least one 
singularity for which \z\ < 1. This singularity implies 
a smaller radius of convergence for the series (|3.53j) . In 
order to check if the function u(z) is convergent at z = 1 
we must analytically continue the series by constructing 
the expansions of u{z) iteratively at some midpoints 67 . 

Namely, we equate the series expansion 



l (z) = j2b n z n = j2~b n (z-z r, 



(3.65) 



n=0 



n=0 



where z = zq is a midpoint inside the radius of conver- 
gence of (|3.53[) . 

The coefficients b n satisfy the ^V-term recurrence rela- 
tion, which could be reduced to the three-term relation 



4? bi + ci} bi-i + c0 6,-2 =0, i > 1. (3.66) 

In order to find 61/62, we must use the condition at the 
event horizon by taking account of (|3.65p . 



s(3)r 



til the cosmological horizon (or spatial infinity) appears 
to be inside the radius of convergence. 

If the convergence of the continued fraction on the 
right-hand side of Eq. (|3.68l) is slow, one can use the 
Nollert improvement. Since the radius of convergence of 
the Frobenius series is now less than 1 (R < 1), we must 



choose Co 



R- 



(l3"62l) . 



K. Generalization of the Frobenius series 



The series expansion (|3.53|) does not need to be neces- 
sarily some power of a rational function of r. For some 
cases, a more convenient choice is an expansion in terms 
of another full set of functions in the appropriate Hilbert 
space. Here we considered, as an illustrative and quite 
general example, charged scalar field perturbations of 
Kerr-Newman-de Sitter black holes that are described 
by the line element 



ds z 



dP_ 



a7 



(l + a)V 



A e sin 2 6 
(1 + a)V 
(dt - a sin 2 9 dtp) 2 



[ctdt- (r 2 + a 2 )dipf 
(3.69) 



60 = V 6„Zq, 61 = V nb n ZQ 



rc=0 



(3.67) where 



We find the coefficients b n from Eqs. (13.58)) and substitute 
them into Eq. p.66[) . If z — 1 is the closest singular point 
to z — zq, we obtain the equation with respect to uj as 



61 
60 



~(3) ~(3)~(3) -(3) -(3) 



~(3)_ -(3)_ -(3) 
c l,4" 



(3.68) 



Otherwise one has to repeat the procedure, constructing 
the series (|3.65[) for the next midpoints zi, 22, Z3, . . ., un- 



A r = (r 2 + a 2 ) (l - -^r 2 ) - 2Mr + Q 2 



1 + a cos 2 1 



Aa 2 



(3.70) 



M is the black hole mass, Q is the charge, a is the rotation 
parameter, and A is the cosmological constant. 

After the separation of variables, the angular part of 
the equation of motion for the massless charged scalar 
field (I2.8[) can be reduced to the following form 2 ^: 



^-(1 + m; 2 )(l - x 2 )4- + A - s(l - a) + H±^!^ _ 2ax 2 + 
dx dx a 

1 + a ( 2s(am - (1 + a)fla - (1 + ^\ 2 - 2m(l + a)g + .s 2 ) - (3.71) 



1 + ax 2 



(l + a) 2 m 2 (1 + a)(s 2 + 2smx)\ 

~ {1 + ax 2 )(l - x 2 ) l^ 2 J ^ 



0, 



where A is the separation constant, £ = auj, x = cosO, s is the field spin, and m is the projection of the angular 
momentum of the field onto the axis of the black hole rotation. Here s and m are (half) integers and < s < 2. 
The appropriate series for the function S is 2 ^ 

00 

S(z) = z Al (z - l) A - (z - z s ) As ( z - Zoo )J2 Ku n (z), (3.72) 

n=0 
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where 



* a. — i x 



i(l + iy/a) 2 
k-sfa. 



A l = 



\m — s\ 



A, 



| m + s | 



2 ' 2 ' 

The expansion is done in terms of the Jacobi polynomials 



i ( 1 + a _ ^„ 
A 3 = ±- (— ^-^-is 



u n (z) = F(-n, n + u,; T ,z) = (-1) ^ + ^ P (fi-^) (2 z 1), 

where Q = 2A X + 2A 2 + 1 and 7 = 2A X + 1. 

The coefficients b n in (|3.72p satisfy the three-term recurrence relation (|3.58al and I3.58b[) with 



„(3) 



„(3) 
-l,n 



.(3) _ 



n(n + Ai + A 2 t s)(n + 2A 2 ) 



2(2n + 2^! + 2A 2 + l)(n + A x + A 2 ) ' 



±^7 



Jn 



+ 



(n - l)(n + 2A 1 + 2A 2 ) 



2(n + Ai + A 2 )(n + Ai + A 2 - 1) ' -I 
A - 2A X A 2 -Ax-A 2 + 2(m + s T (2^i + 1))? - 



1 

i (71 - 1 + Ai + A 2 ± s)(n - 1 + 2Ai)(n - 1 + 2A X + 2A 2 ) 
T v^ € 2(2n + 2A 1 + 2A 2 -3)(n-l + A 1 +A 2 ) 



(3.73) 
(3.74) 

(3.75) 



where 



J n = {n- l)(n + 2A X + 2A a )(Ai - A 2 ) + {A x + A 2 ± s + 1) ( (n - l)(n + 2^ + 2A 2 ) + (2A X + + A 2 ) ) . 

I 



Since the series Q3.72p must be convergent at z = 1, 
we can solve Eq. (|3.60[) numerically in order to find the 
separation constant as a function of frequency. 

In the Reissner-Nordstrom-de Sitter limit (a — > 0) 

Eq. (|3.60j) is reduced to cp„ = 0. In this case the value 
of A does not depend on to 

A = (£-s+l)(l+s), t = n+Ai+A 2 > max(\m\, \s\). 



L. Example I: Frobenius series for radial part of the charged 
scalar field in the Kerr-Newman-de Sitter background 

The radial part of the massless (charged) field equation 
of motion is of the form (|3 . 5 1 1) and reads22 



The appropriate Frobenius series is found to be 



R(r) 



r+ 



-s-2iK(r+)/A' r (r+) 



e: B(r) r -2s-l u 



• + 1 00 



r - r_ Too - r + 



(3.76) 



2a 



(s + l)(2s + l)r 2 + 2s(l - a)-2isqQ~ A } R(r) = 0, 



Note that in order to obtain the recurrence relation for 
both types of the boundaries (asymptotically flat and 
de Sitter), we introduced the exponent e* B ^ such that 
dB(r) K mi . 

— - — = — — . This exponent describes outgoing waves 
ar A r 

at horizons and spatial infinity. Thus, we have to com- 
pensate the outgoing wave at the event horizon. That is 
why the factor of 2 appears in the power of the first mul- 
tiplier in p.76[) . The parameter ro is fixed to be the inner 
horizon r_ in order to move all the singularities outside 
the unit circle and, at the same time, to provide the best 
convergence of the infinite continued fraction (|3.60l) . 

In the Reissner-Nordstrom-de Sitter limit (a = 0) for 
the uncharged field (q = 0) we obtain K = cor 2 and 



where K — [uj(r 2 + a 2 ) — am](l + a) — qQr, and q is the 
field charge. 



R(r) 



— s—iuj/n 

r ~ r + \ iuir t -2s- 1 

e r I 



r - r + ^ - r_ 



r — r_ r„ 



(3.77) 
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The tortoise coordinate is denned as dr* = 



r dr 



2r\ 



A, 



and 



is the surface gravity on the event horizon. 



If A > 0, we can observe the same asymptotic behavior 
of the exponent at the cosmological horizon as for the 
case of fj, — 0: 

iA(r) 



M. Example II: Frobenius series for the massive scalar field 
in the higher-dimensional Reissner-Nordstrom-de Sitter 
background 

The D-dimensional Reissner-Nordstrom-de Sitter 
black hole is described by the metric 



dr 2 

ds 2 = -f(r)dt 2 + T -^+ r 2 dVt D 

fin 



(3.78) 



where cK7£>_ 2 is a metric of a (D — 2)-dimensional sphere, 
2M Q 2 2Ar 2 



fir) = 1 - 



(D-1)(D -2) 



After separation of the angular and time variables the 
radial part of the massive scalar field equation of motion 
(□ — n 2 )^ = is reduced to the wavelike equation 



— + u 2 - f(r)V{r) ) r^R(r) = 0, (3.79) 



where the effective potential is 



V{r) = ^ 



2 M+D-3) J'(r)(D-2) J(r)(D-2)(D-4) 



2r 



4r 2 



where i parametrizes the angular separation constant. 
The Frobenius series for this case is 



r ) __-__ (r -r + roc - R 



r - R ^ - r + 



RM = ( ____± \ e *Mr) r , „ 
r — R 

where k = ~/'(r + ), e lA ^ describes the outgoing wave 
for spatial infinity and the horizons and satisfies 



dA(r) y/uj 2 - fi 2 f{r) 



dr 



fir) 



The sign in the exponent is fixed by the quasinormal 
boundary condition: The real part of A(r — > oo) must 
be of the same sign as the real part of the eigenfrequency 
bj. This choice of the sign makes the wave outgoing at 
spatial infinity. 

For the massless field (/_ = 0) this exponent is e" 1 ""*— 
and the Frobenius series reads 



,r-r + \ K 4 _e-2 (r -r+ r^ - R 
Rir) = I ir I e *r 2 u 



r-R 



r - R Too - r i 



Since for D — Awe can choose R = r_ , we come to (|3.77[) 

(8 = 0). 



For the asymptotically flat case the asymptotic behav- 
ior is different in four- and in higher-dimensional space- 
times. If D > 5, f(r) = 1 + o(r _1 ) and the Frobenius 
series can be determined as£& 



Rir) 



r-R 



,>„ A .2-i/2 p-2 r — r + 
e' r V" ^ r ~u 1 



r-R 



For D = 4 the term of order ~ 
nontrivial contributio n 48 ' 49 



(3.80) 

in f(r) leads to the 



R(r) = 



R 



(3.81) 



r 2tM^^ +l M^/y^JJ r -^ u ( r-r+ 

r — R 



The same approach could be applied for the Maxwell 
field and the gravitational perturbations, because the ra- 
dial parts of their wave equations can be reduced to the 
form d33Ip3-^. 



N. Horowitz-Hubeny method 

In order to find quasinormal modes in the asymptoti- 
cally anti-de Sitter space-times we normally need to im- 
pose Dirichlet boundary conditions at the spatial infinity. 
Thus, we find the appropriate expansion for the function 
R(r) in (|3.51[) without consideration of the singularity 
point at the infinity. This method was proposed ir>2_.. 
Namely, we define 



R = z-^l (2 ^i>(z), 



(3.82) 



where k is the surface gravity at the event horizon, z — 
r — r~L- _____ _____ 

. Substituting (I3.82[) into (|3.51[) one can rewrite 

r — ?*_ 

the radial equation in the form 

a(zW(z) + -^V(z) + = 0, (3.83) 

z z A 



s(z) = ^s n z", t(z) = ^2t n z n , u(z) = ^2u n z n , 



with respect to the wave function ip(z) = a n z n , which 

is regular at the event horizon 2 = 0. The Dirichlet 
boundary condition ip(z = 1) = implies 



(3.84) 



n=0 
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The coefficients a n can be found through the recurrence 
relation 

Sfe=o a k(k(k - l)s n -k + kt n -k + Un-k) fo oc n 
a n — -, 7T ; — t 7 (3.8oJ 



starting from an arbitrary ao- Substituting (|3.85p into 
Eq. (|3. 841) . we find the equation with respect to eigen- 
value uj. Since the sum (|3.84[) is convergent, one can 
truncate the summation at some large n and minimize 
its absolute value with respect to uj. Doing so we in- 
crease the number of terms in the sum until the value of 
uj does not change within required precision. 

Note that for Eq. (|3.84p to be convergent, r_ has to 
be a singular point, and all the other singularities, ex- 



cept r = r + and r = 



must lie outside the unit circle 



\z\ > 1. If both of these conditions are impossible to 
satisfy simultaneously, we must use the continued frac- 
tion method with the appropriate fixing of the behavior 
of R at spatial infinity. The Horowitz-Hubeny method 
converges slower than the Frobenius one. 



O. Limit of high damping: Monodromy method 

The limits of high damping for QN spectra of asymp- 
totically flat or de Sitter and asymptotically anti-de Sit- 
ter black holes have qualitatively different behavior. For 
asymptotically flat or the de Sitter case 

|Im(w)| > |Re(w)|, n^oo (3.86) 

while for asymptotically AdS black holes 

|Im(w)| ~ |Re(w)|, n ->• oo. (3.87) 

The monodromy technique starts by considering the wave 
equation (|2.12[) in the complex plane, so that instead of 
the real tortoise variable r* one uses the complex variable 
x. Then, Eq. (12.12)) can be written in the form 



d 2 *(r) 
d(ujx) 2 



+ (1 - K(r)w _2 )#(r) = 0. 



(3.88) 



In the limit \uj\ — > oo the term proportional to uj~ 2 in 
Eq. (|3.88p can be discarded, so that the general solution 
can be written as a superposition of plane waves 



= Ae +lux + Be 



(3.89) 



Then one needs to find the Stokes line as contours in the 
complex ujx plane. The Stokes lines are defined as lines 
satisfying the equation lm(ujx) = 0. If we apply Eq. (I3.4[) 
and x = XRe — i%im-> then 

LUX = UJ Re X Re - UJ Im Xl m - i(uJ Im X Re + UJ Re Xl m ) = 

UR e x Re - uj Irn x Im , lm(ujx) = (3.90) 

along the Stokes lines. Thus, along the Stokes line, the 
plane wave solution (|3. 891) oscillates without decay (or 



growth). The next step is to use the basic numerical 
data for quasinormal modes with large overtones, such 
as Eq. p. 861) or (|3.87p in order to translate the condi- 
tion for the Stokes lines in the ojx-plane to the condition 
in the complex x-plane. For instance, for asymptotically 
flat or de Sitter black holes the inequality f|3 . 86[) imme- 
diately gives Re(u;) = 0, which is the anti-Stokes line in 
the complex ir-plane. Then one needs to determine the 
closed contour in the complex x-plane and calculate the 
monodromy, which will give an analytic expression for 
the quasinormal frequency uj. Although the approach is 
well known from complex analysis, in the present form 
it was used for the problem of highly damped quasinor- 
mal modes only recentl y 69 i 70 for asymptotically flat cases 
and developed ir*Zir— for various asymptotically dS and 
AdS black holes. Particular elements of the monodromy 
method, the Stokes lines, were used for QNMs ir4£. 



IV. QUASINORMAL MODES: SUMMARY OF RESULTS 

A. More on the context in which quasinormal modes are 
studied 

The quasinormal spectrum of a stable black hole is 
an infinite set of complex frequencies which describes 
damped oscillations of the amplitude^. It is clear that 
the oscillation with the smallest damping rate is domi- 
nant at late time, while oscillations with higher damping 
rate are exponentially suppressed. In practice, by study- 
ing the signal at the stage when quasinormal ringing is 
observed, we are able to extract a few dominant modes 
which have the smallest imaginary parts. 

The quasinormal ringing was first described by 
Vishveshwara when considering the scattered packet 
of gravitational waves in the background of the 
Schwarzschild black hole^. The term quasinormal fre- 
quencies was introduced by Press in 19712S. In the 
same year the lower quasinormal modes were calculated 
by studying test particles falling into the Schwarzschild 
black holeZ£. Later, quasinormal modes were calculated 
for Kerr black holes--. The quasinormal spectra of the 
coupled gravitational and electromagnetic perturbations 
of the Reissner-Nordstrom black holes were first studied 
inZ&. 

From that time until now quasinormal modes have 
been extensively studied within various contexts, as fol- 
lows: 

1. Ringing of astrophysical black holes: In this 
context the dominant frequencies of the quasinor- 
mal spectrum were studied for astrophysically mo- 
tivated black holes in four dimensional asymptoti- 
cally flat and de Sitter space-times (see Sec. IIV.C1) . 

2. Analysis of stability and quasinormal modes 
of mini black holes: This is an important issue in 
various extra dimensional scenarios, such as ADD 80 
and the Randall-Sundrum 81 models. In this con- 
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text dominant quasinormal frequencies of higher- 
dimensional black holes, black strings and br aires, 
Kaluza-Klein black holes and other black objects 
have been extensively studied (see Sec. IIV.D[) . 

3. In the context of the AdS/CFT correspon- 
dence: Here the quasinormal modes of large AdS 
black holes and of £>p-branes are studied (see 
Sec. El. 

4. Loop quantum gravity interpretation: It was 

suggested that the asymptotic real part of quasinor- 
mal modes of black holes can be interpreted within 
loop quantum gravity, which allows us to fix the 
Barbcro-Immirzi parameter, the free parameter of 
the theory (see Sec. IVf . 

5. Lower-dimensional black holes: These were 
studied mainly because of their simplicity, so that 
in many cases it is possible to find exact solutions 
for the quasinormal spectra. For instance, it was 
shown that the quasinormal modes of the (2 + 1)- 
dimensional AdS black hole exactly coincide with 
the poles of the retarded Green's function for the 
two-dimensional CFT (see Sec. IIX.BI) . 



B. Isospectrality 



Consider two wavelike equations (|2.12j) with the effec- 
tive potentials V + and V~ such that 

dW{n) 



V ± {r i ,) = W 2 (r*)± 



dr+ 



(4.1) 



where W(r+) is some finite function, ft is a constant, and 
r* is the tortoise coordinate. If ^ + is an eigenfunction 
of the wavelike equation with the potential V + , then the 
eigenfunction for the potential V~ is given by 



*"(r*) cx 1%) - - * + (r*), (4.2) 

which corresponds to the same eigenvalue u). Therefore, 
the quasinormal spectrum is the same for the potentials 
V+ and 



If one takes 
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then V + and V~ are the effective potentials for gravita- 
tional perturbations of axial and polar types respectively 
for the Schwarzschild black hole2£, 



V+ = f(r] 



f £(£+1) 6M\ 



r 3 
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FIG. 5 The first 60 quasinormal modes for the gravitational 
perturbations of the Schwarzschild black hole; the numerical 
values of 1000 lower quasinormal modes are available from 



http://qnms.way.to 



This means that the quasinormal spectra of gravitational 
perturbations of axial and polar types coincide. The 
same symmetry is preserved in a more general case of 
the Reissner-Nordstrom-de Sitter black hole, where both 
types of gravitational perturbations are coupled to elec- 
tromagnetic perturbations. If, for instance, one considers 
a scalar field which is coupled only to the polar pertur- 
bations, the isospectrality is broken. 

The isospectrality between axial and polar perturba- 
tions does not exist for higher than four dimensional 
black holes, when, in addition to the latter two types of 
perturbations, the other dynamical type appears. This 
new type of gravitational perturbation transforms simi- 
lar to a tensor with respect to coordinate transformations 
on a (D — 2)-sphere. It is interesting to note that the 
tensor-type gravitational perturbations are described by 
the same effective potential as a test scalar field in the 
black hole background. 



C. Quasinormal modes of D — 4 black holes 

The quasinormal spectra of four-dimensional black 
holes were studied starting from the 1970s in the con- 
text of possible observations of ringing from astrophysi- 
cal black holes (see, e. g-r&). Figure [S] shows the lowest 
gravitational quasinormal modes of Schwarzschild black 
holes, calculated accurately with the help of the contin- 
ued fraction method (see Sec. lIII.G|) . The imaginary part 
of the higher overtones grows, while the real part de- 
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creases until it becomes almost zero at some moderately 
large overtone number n (for £ — 2, n — 9, ior £ = 3, 
n = 41). The corresponding purely damped mode is ap- 
proximately equal to the so-called "algebraically special 
mode" , which is given by§^ 

uM = -i{£ -!)£{£ +!){£ + 2) /12. (4.3) 

For modes which are higher than the algebraically spe- 
cial one, the real part starts growing and approaches 
its asymptotic value, which can be found exactly (see 
Sec. IIV.F|) . Since the damping rate of the algebraically 
special mode (|4.3[) grows quickly with £, the asymptotic 
regime is achieved at high overtones. For instance, for 
I — 6 the asymptotic regime of high damping is achieved 
atn - 10 5 S. 

Real astrophysical black holes are rotating. In addi- 
tion, one can neglect their electric charge so that the 
Kerr metric is the most astrophysically motivated exact 
solution of the Einstein equations. Thus, gravitational 
quasinormal modes of Kerr black holes are of primary in- 
terest for observations of gravitational waves. Accurate 
calculations of quasinormal modes for Kerr black holes 
were performed by Leaver—. The basic properties of the 
spectrum for small rotation were found from the eikonal 
limit within the slow-rotation approximation. Using the 
Poschl- Teller formula ([3~T2]) . one can find thal£i 

t » \m\ > 1, a < M, (4.4) 

where m is the azimuthal number and a is the rotation 
parameter of the black hole. 

Rotation gives us the following new properties of the 
quasinormal spectrum: 

• The real part of the quasinormal mode grows with 
am. 

• The quasinormal modes explicitly depend on the 
azimuthal number m and there is a symmetry 
uj(rn) —¥ — w*(— m). This allows us to consider only 
modes with the positive real part. 

Table U summarized the publications in which the 
quasinormal ringing of four-dimensional black holes was 
studied with the help of various numerical methods. The 
reader is referred to the original works (see Table [XJ) for 
more details and numerical data. In Sec. IIV.HI we dis- 
cussed dependence of the quasinormal modes on the black 
hole parameters. 

Lower modes of the Dirac field were found by the 
WKB approach for SchwarzschilcU^ and Schwarzschild- 
de Sitter backgrounds 1 ^ and by using the approximation 
with the Poschl- Teller potential (see Sec. IIII.BI) for Kerr- 
Newman-de Sitter black holes 104 . Later the spectrum of 
the Dirac field was calculated using the accurate Frobe- 
nius method for Schwarzschild 105 , Reissner-Nordstrom- 
de Sitter^S 6 -, Kerr-Newmar>i2£ and Kerr-Newman-de Sit- 
ter black holes 1 ^. The quasinormal spectrum of the 



TABLE I Publications where the lower modes of various four- 
dimensional black holes were calculated by different methods: 
integration of the equations of motion, the WKB formula, and 
the method of continued fraction. 



QNMs 


Integration 


WKB 


Cont. f. 


Schwarzschild 


77 


85 


64 


Reissner-Nordstrom 


79 


86 


87 


Reissner-Nordstrom (extreme) 






88 


Kerr 


78.89 


90 


64 


Kerr-Newman 




91 


55 


Reissner-Nordstrom+dilaton 




92.93 




Topological black holes 


94 






Stringy black holes 




95 




Schwarzschild-de Sitter 


96 


97 


98 


Reissner-Nordstrom-de Sitter 


99 


99 




Kerr-de Sitter 






100 


Garfinkle-Horowitz-Strominger 




101 





scalar and Dirac fields around the Born-Infcld black hole 
was studied m 1 ^. Dominant gravitational frequencies of 
black holes in the scalar-tensor gravity were calculated 
mil . 

Real astrophysical black holes are not isolated, but 
surrounded by some matter. Theoreticians frequently 
call such black holes with an "environment" as "dirty" 
black holes. These were studied using a perturbative 
method iniii. Higher modes of the spherically symmetric 
dirty black holes were studied ir>ii2. Quasinormal modes 
of black holes surrounded by quintessence were studied 
1^113-^116 ; anc [ moc j es f t ne phantom scalar field ioiil. 

Recently, perturbations of more exotic objects such 
as wormholes, white holes, or naked singularities have 
been investigated. It was found that the quasinor- 
mal spectrum of Schwarzschild black holes differs from 
wormholes 1 ^ and white holes 1 ^. Therefore, these ob- 
jects, if they exist, might be detected through observa- 
tions of their quasinormal ringing. 

D. Quasinormal modes of mini black holes 

Quasinormal spectra of black holes attracted consider- 
able interest in the following extra dimensional models: 

1. The large extra dimensions scenario allows for 
the size of extra dimensions to be of a macroscopic 
order—. When the size of the black hole is much 
smaller than the size of the extra dimensions, the 
black hole can be considered as effectively living in 
a D-dimensional world and, thereby, approximated 
by a solution of higher-dimensional Einstein equa- 
tions. The simplest example of such a solution is 
the Tangherlini metric^ 6 ., which is a generalization 
of the Schwarzschild metric for D > 4. 

2. Randall-Sundrum models* 1 assumed that our 
world is a brane in higher-dimensional anti-de Sit- 
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ter warped space-time. The AdS space-time implies 
quick decay of the fields outside the brane. The 
warp factor, which is the parameter of the theory, 
can be set up in order to obtain large size (of the 
order of TeV) or, if one wishes, small size of the 
extra dimensions. In the latter case the gravita- 
tional perturbations contain Kaluza-Klein modes, 
which share many properties of massive fields. For 
instance, the quasiresonances were observed in the 
spectrum of black strings 1 ^ -. The ringing of the 
Randall-Sundrum braneworld was studied i n 121 i 122 , 
where quasibound states of the gravitational per- 
turbations between the AdS boundary and the 
brane were observed. 

3. Brane-localized fields: When the size of extra 
dimensions is larger than the size of the black hole, 
one can consider the model, for which the standard 
model particles are restricted to live on a (3 + 1)- 
brane, while gravitons propagate also in the bulk. 
When considering evolution of the standard model 
fields in the background of a mini black hole, one 
can think that the mini black hole effectively be- 
haves similar to a higher-dimensional one projected 
onto the brane 1 ^ 3 -. 

Quasinormal modes of a scalar field in the background 
of a D-dimensional black hole were estimated ini 2 ^ using 
the third-order WKB approach and more accurate values 
were obtained inii with the sixth-order WKB formula 
(sec Scc. IIII.Dj) . The scalar field in the background of the 
five-dimensional rotating black hole was studied in 1 ^. 

Quasinormal modes of gravitational perturbations of 
the Tangherlini black hole were calculated in 42 using the 
higher-order WKB formula. Later, five-dimensional 126 
and higher-dimensional 6 Tangherlini black hole pertur- 
bations were studied using the Frobenius method. Fi- 
nally, quasinormal modes of all types of perturbations for 
higher-dimensional Reissner-Nordstrom-de Sitter black 
holes were calculated i n 68 ' 127 . Quasinormal modes of the 
Dirac field in the background of higher-dimensional black 
holes were found ir>i2£. 

The quasinormal spectrum of a scalar field in the back- 
ground of ultraspinning black holes was studied in six-J^ 
and higher-dimensional space-times^ -. In D > 7, quasi- 
normal modes of tensor type of gravitational perturba- 
tions were studied for simply rotating Myers- Perry black 
holesi 3 - 1 . and their asymptotically AdS generalizations 1 ^. 
It was found that the rotating AdS black holes suffer from 
the superradiant instability (see Sec. IVIII.F[) . 

For mini black holes the quantum corrections are 
significant 1 ^. These corrections are suggested by the 
string theory. Its slope expansion yields higher-order cur- 
vature corrections to the Einstein action, which make the 
gravity be mathematically non-contradictory in higher 
dimensions. The first order correction to Einstein grav- 
ity is proportional to the square of the curvature and 
is a topological invariant in four-dimensional space- 
time 1 ^. This correction, called the Gauss-Bonnet term, 



implies a nontrivial correction to the higher-dimensional 
backgrounds 1 ^. The quasinormal spectrum of the scalar 
field was studied for the neutral^ 6 - and charged 1 ^ 7 . Gauss- 
Bonnet black holes, as well as for their asymptotically de 
Sitter and anti de Sitter generalizations 1 —. 

The quasinormal modes of tensor- and vector-type 
gravitational perturbations of Gauss-Bonnet black holes 
were calculated using the third-order WKB formula. 
The quasinormal modes of gravitational perturbations of 
all types and the development of an instability in time 
domain were found in 1 ^ -. 

The ringing of the brane-localized standard model 
fields was studied in the background of Schwarzschild, 
Schwarzschild-(anti) de Sitter^ 41 -, Ker r 142 ' 143 , and Gauss- 
Bonnet black holesi 4 ^. The standard model fields were 
also considered in the scenario of the three-brane with fi- 
nite tension in the six-dimensional bul k 145 i 146 and on the 
two-brane in the four-dimensional version of the Randall- 
Sundrum modet^. 

In higher-dimensional space-time one can consider 
more complicated black hole configurations, which were 
studied for the particular case of five-dimensional space- 
time. When a gauge field is added, a five-dimensional 
black hole may be embedded in a rotating universe, which 
is called the Godel universe. The quasinormal modes 
of the scalar field in the background of such a black 
hole in the Godel universe were studied ini 4 ^. Another 
family of solutions in the five-dimensional space-time is 
the squashed Kaluza-Klein black holes, which asymptot- 
ically look similar to four-dimensional black holes with 
the Kaluza-Klein modes due to the compactified extra di- 
mension. The quasinormal spectrum of such black holes 
allows us to probe the extra dimensions even at small 
energies, for which the Kaluza-Klein modes are not yet 
excited. The quasinormal modes of scalar field and grav- 
itational perturbations of squashed Kaluza-Klein black 
holes were calculated ini 4 ^ and confirmed by indepen- 
dent calculations oS^°-. Then, an analysis was continued 
to the squashed Kaluza-Klein black holes in the Godel 
universe 1 ^. 



E. Quasinormal modes of AdS black holes 

In asymptotically AdS backgrounds there is no outgo- 
ing wave at spatial infinity. The AdS boundary provides a 
confinement for perturbations, which are, thereby, local- 
ized relatively closely to the black hole. The quasinormal 
ringing governs evolution of perturbations at all times, 
i. e., there is no late-time tail stag o 30 ' 155 as in the asymp- 
totically flat and de Sitter space-time (see Sec. IVII|) . The 
high overtone regime is usually already reached at mod- 
erate overtones n, and one can find that the real and 
imaginary parts of asymptotically high overtones have 
equidistant spacing. 

Within the AdS/CFT correspondence the quasinormal 
modes of black holes in the anti-de Sitter space-time have 
a direct interpretation in the dual conformal field the- 
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ory (see Sec. IIXI) . The quasinormal modes of the con- 
formally invariant scalar field in the Schwarzschild-AdS 
background were first studied ii*^ prior to the inter- 
est in the ADS/CFT correspondence. In the context of 
the AdS/CFT correspondence, quasinormal modes of a 
test scalar field were calculated in^. Later, quasinor- 
mal modes of the electromagnetic field and the grav- 
itational perturbations of the Schwarzschild-AdS black 
hole were found fr>i£2. Quasinormal modes of toroidal, 
cylindrical and planar black holes in anti-de Sitter 
space-times were calculated ir*i£l. Quasinormal ringing 
of the Reissner-Nordstrom-AdS background was stud- 
ied for scalar fielcU^, gravitational and electromagnetic 
perturbation*^ 6 -. 

When the radius of the black hole is much smaller than 
the AdS radius, the effect of the black hole is small and 
the quasinormal modes approach the normal modes of 
the pure AdS space-tim o 157 i 158 i 223 . For large black holes 
the real and imaginary parts of quasinormal modes grow 
proportionally to the black hole radius. The only excep- 
tion is the fundamental vector and scalar mode o 42 i 156 i 159 , 
which are hydrodynamic modes (sec Sec. ILX.Cj) 160 . 

Quasinormal modes of the Dirac field were found 
for the Schwarzschild-AdS 161 and Reissner-Nordstrom- 
AdSi^ 2 - black holes. In addition, quasinormal ringing 
was studied for the four-dimensional Kerr-AdS black 
hole^ 3 - and planar AdS black holes in four and higher 
dimensionsi^-i 6 ^. An alternative, perturbative ap- 
proach for finding quasinormal modes was developed 



Schwarzschild-de Sitter black hole 
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,167.168 



for asymptotically flat and AdS black holes. 



Using an analytic technique based on the complex co- 
ordinate WKB method, the so-called "highly real" quasi- 
normal modes were reported to appear in the spectrum of 
Schwarzschild-anti de Sitter black holes. The real part of 
these quasinormal frequencies asymptotically approached 
infinity while the damping rate remained finit o 169 ! 170 . 
Since for large black holes the imaginary part of highly 
real quasinormal modes is less than those which are 
established as the fundamental modes, the highly real 
quasinormal modes must be dominant in the spectrum. 
Nevertheless, these highly real modes were not observed 
through time-domain integrations nor any alternative 
calculations. 



F. Exact solutions and analytical expressions 

Here we summarize some analytical formulas for quasi- 
normal modes. 

• Large multipole number: Quasinormal modes 
can be obtained by substituting the corresponding 
effective potentials into the first-order WKB for- 
mula and taking the limit of large £. 

The large multipole limit for the Schwarzschild 
black hole was derived in 3 - 7 - and generalized for the 
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1 



_ VI - 9AA/ 2 

~ :'.\7i.U V -* V" ' 2 
where n is the overtone number 



O 



(4.5) 



This formula is valid also for scalar (£ = 0, 1, 2, . . .), 
electromagnetic (£ = 1,2,3,...) and Dirac (£ — 
1/2, 3/2, 5/2,...) fields. Despite the fact we consid- 
ered large £ expansions, the formula (I4.5[) also gives 
good estimations for the dominant frequencies with 
small £. For instance, for the £ = 3, n = gravita- 
tional perturbation of the Schwarzschild black hole 
(A = 0), one has u> = (0.599 - 0.093i)M -1 , while 
formula (JO} gives U) = (0.674 - 0.096i)M _1 . 

For higher-dimensional black holes the large multi- 
pole expansion reads^ 



1 



4 

21 + D - 3 - 



D 



D - 1 



(4.6) 



2n + 1 



This was generalized to the multihorizon black 
holes 
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• Asymptotically high overtones: For 

Schwarzschild black holes, asymptotically high 
overtones satisfy the relation 6 - 9 - 



e 8 ™ M = -l-2cos(irs), 



Im(w) < 0, 



(4.7) 



where s = 0,1,2 for the test scalar, Maxwell, and 
gravitational fields, respectively. The solution to 
Eq. (gUJ) is 



2Mw 



ln(3) n-l/2 
Itt 2 



2 ' 



i, scalar, gravitational; 



Maxwell. 



(4.8) 

The asymptotic behavior of quasinormal modes of 
the Dirac field is the same as for the Maxwell field: 
the real part approaches zero, while the spacing of 
the imaginary part tends to 1/ 4M-ifl&. 

Formula (|4.7[) was generalized for arbitrary single- 
horizon black holes asi 7 - 2 - 



T — 2cos(7rs) 



Im(w) < 0, 



(4.9) 



where Tjj is the Hawking temperature. 

The analog of formula (|4.7j) was derived also for 
charged black holes, for black holes in the asymp- 
totically de Sitter space-time, and for higher- 
dimensional black holes (see Table HI)) . Yet, these 
equations cannot always be solved algebraically 
with respect to w (see topical review-Z 3 -). 
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TABLE II Publications where the higher modes of various 
black holes were found. 



Black hole 


Asymptotic formula 


S chwarz sch ild 


^ (numerical),^ 


Dirac field in Schwarzschild 


105 




Tangherlini 


70,173 174 






(numerical) 


R,eissner-Nordstrom 


33. {D = 


4),B (D > 4) 


Kerr 


iffi (numerical) aLill 


J. vCI 1 11CW llloJ.1 


178 




Spa 1 fit fipld in 1VT vpts- Pprr v 


179 




Single- horizon (e. g. , dilatonic) 


172 




Spherically symmetric 


aa (inte 


ger spin),— 


Garfinkle-Horowitz-Strominger 


iSa (scalar),^ (Dirac) 


Gibbons- IVlaeda dilaton 


184 




Gauss-Bonnet 


185 




i, n pt*T1 — Sim Pin d 

V_>±±C1 11 VJllllUJ-lu 


186 




Schwarzschild-de Sitter 


— (numerical)^ 


TVin P'hprlini-H p Siffpv 

_1_ CLlltilld ±1111 VIC kJl Ij Lj CI 


73 




R,eissner-Nordstrom-( A)dS 


^ (numerical),- 2 (-D > 4) 


EM field in SdS (D > 4) 


188 




Extremally charged RN-dS 


189 




Large SAdS (D = 4) 


72 




Large SAdS (D = 5) 


^o. (numerical) ^SLiSS. 


All SAdS 


73 




Massive scalar field 


48 (D = 


(D > 4) 



• Near extreme black holes in the de Sitter 
space-time: When the black hole horizon ap- 
proaches the cosmological horizon, the effective po- 
tential approaches the Poschl- Teller potential (see 
Sec. IIII.B[) and the quasinormal spectrum can be 
found analyticall y 35 ' 36 . For the four-dimensional 
Schwarzschild-de Sitter black hole the spectrum has 
the form 

~ = ^+i- s )(£ + s )-I - (n + ~) i, (4.10) 

where s = 0, 1,2 for the test scalar, Maxwell, and 
gravitational fields, respectively, £ = s, s + 1, . . . is 
the multipole number, and k is the surface gravity 
at the event horizon. The above formula was later 
confirmed with the help of the accurate Frobenius 
method^. 

• Normal modes in the pure anti de Sitter 
space-time (without a black hole): The effective 
potential in this case is an infinite potential well, 
where bound states exist. There is nothing that 
could absorb the energy of these states. This is why 
the quasinormal spectrum consists of normal modes 
only, i. e., the imaginary parts of the frequencies are 
zero. The quasinormal spectrum was first found 
ini2£ for a massless scalar field. For gravitational 
perturbations of four-dimensional AdS space-time 



the spectrum was calculated ini 5 ^ and generalized 
for an arbitrary number of space-time dimensions 

• 73 

as 

LdR=2n + D + £-j, n G N, 

where R is the anti-de Sitter radius, D is the num- 
ber of space-time dimensions, I is the multipole 
number, j = 1 for the scalar field and the grav- 
itational perturbations of tensor type, j — 2 and 
j = 3 for vector and scalar types of gravitational 
perturbations. 

• Quasinormal modes in the de Sitter universe: 

Perturbations in a de Sitter universe decay as en- 
ergy can pass through the cosmological horizon. 
The quasinormal spectrum differs in spaces with 
odd and even D and can be described by different 
formulas. The exact solution for odd D was found 
in^ 3 - and does not depend on D 

UToo = -i(2n + £ + j - 1), n <E N. 

Here, r m is the cosmological horizon. For even D 
the spectrum has two branches^ 5 - 

ujtoo = — «(2n + I + j - 1), neN 

and 

ojtoo = -i(2n + D + £- j), neN. 

It is interesting to note that the second branch 
of quasinormal modes coincides with the normal 
modes of AdS space-time if we substitute an imag- 
inary value of the radius of the cosmological hori- 
zon. 

• Quasinormal modes of the Schwarzschild 
and Kerr black holes: The recent development 
of the theory of Heun equations allows us to find 
exact solutions for the wave equations in four- 
dimensional space-time. The Regge- Wheeler equa- 
tion (|2.13[) . which describes perturbations of the 
Schwarzschild black hole, can be reduced to the 
confluent Hcun equation and the quasinormal spec- 
trum can be expressed as a solution of nonalge- 
braic equations with Heun functions. This equa- 
tion can be solved numerically showing an excellent 
agreement with the known results^ 6 -. The pair of 
equations which describes perturbations of the Kerr 
black hole for any s have been recently solved in-^ 
as well. 

Despite the fact that the quasinormal spectrum for 
such four-dimensional black holes can be expressed 
as a solution of a set of equations, the calculation 
of the quasinormal spectrum remains a nonelemen- 
tary procedure because of dealing with a much less 
studied class of special functions. 
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• BTZ black hole is the (2 + l)-dimensional so- 
lution of the Einstein equations in asymptotically 
anti-de Sitter space-time^ 8 . In 2 + 1 space-time 
there are no gravitational degrees of freedom, al- 
though there are dynamical perturbations of test 
fields. The equations of motion for the test scalar 
and electromagnetic fields are the same and can be 
solved analytically in terms of hypergeometric func- 
tions. The quasinormal spectrum has the formic 

UjR = ±m-2iVM(n + l), (4.11) 

where R is the AdS radius, M is the black hole 
mass, and m is the azimuthal number. 

Since in 2 + 1 gravity the Riemann tensor is com- 
pletely determined by the matter source, the only 
radiative correction to the geometry comes from 
quantum excitations of the matter fields, and the 
perturbative expansion receives no corrections from 
graviton loops. Thus, it is possible to find the 
one-loop quantum correction to the BTZ black 
hole metric and the corresponding correction to the 
quasinormal spectrum, which is given by 200 

_ 2 _ 2 4m 2 + Ml p F(M) 

where lu is the frequency of the BTZ black hole, 
l p is the Planck length, and F{M) is a nonlinear 
function which exponentially approaches zero for 
large M22I. 

For the rotating BTZ black hole the exact expres- 
sion for the quasinormal spectrum of a massive 
scalar field was also found 2 ^ 2 - 

jR = ±m -2i r -±±^ L + 1 + ^+~» y (4.12) 

where /i is the mass of the scalar field, and r_ and 
r + are the inner and outer horizons, respectively. 

Note that for the (2 + l)-dimensional de Sitter ro- 
tating black hole there are several families of the 
quasinormal frequencies, which are given by (I4.12p 
with imaginary values for R and r_ . Some of those 
modes have positive imaginary parts, implying that 
the black hole is unstable^. 

• Other exact solutions: The analytical expres- 
sions for the quasinormal spectra of topological 
AdS black holes were found for the massive scalar 
field 204 and for gravitational perturbations^ - 5 -. 
Also the exact solutions for quasinormal modes of 
dilatonic black holes were found hr 206 ' 207 . 



G. Quasinormal modes of massive fields 

Consider a massive scalar field in the space-time of 
an asymptotically flat black hole. From (|3.80p one can 



see that if one requires a purely outgoing wave at spatial 
infinity, the behavior of the scalar field amplitude at large 
distance behaves as 

where fj, is the mass of the field. 

This implies that for wj m ~ and to < fi there is 
no energy transition to infinity. If, for some values of 
the black hole parameters, the ingoing wave amplitude 
at the event horizon is much less than the amplitude far 
from the black hole, there is no leak of energy from the 
system and one can observe an analog of standing waves. 
An almost purely real mode appears in the quasinormal 
spectrum and the oscillations of the particular frequency 
become long lived. This phenomenon is called quasires- 

49 

The quasinormal spectra of a massive, neutral scalar 
field for Schwarzschild and Kerr black holes^ and for a 
massive charged scalar field in the vicinity of a Reissner- 
Nordstrom black hole 2 ^ were studied first within the 
WKB approximation. It was found that as the field mass 
increases, the damping rate becomes smaller, which was 
confirmed by time-domain integration 210 . Then, accu- 
rate massive modes were obtained in 4 ^ using the conver- 
gent Frobenius method (see Sec. IIII.GI) . It was shown 
that increasing the mass of the field or the mass of the 
black hole gives rise to decreasing of the imaginary part 
of the quasinormal modes until the fundamental mode 
reaches a vanishing damping rate. When some thresh- 
old value of the mass is exceeded, the fundamental mode 
"disappears" from the spectrum and the first overtone 
becomes the fundamental mode (see Fig. [6J . Since this 
mode has a high damping rate, we observed a kind of dis- 
continuity of lifetime values of the oscillations: For some 
particular values of the field mass the oscillations almost 
do not decay, while for slightly higher mass we see that 
the new fundamental mode decays quickly. 

The effect of the field mass is crucial for the lower 
quasinormal modes only. The asymptotic behavior of 
highly-damped modes does not depend on the field 
mass 48 and is the same as for the massless field. 

Quasiresonance for the Kerr black hole was observed 
in^ii, where the stability of the quasinormal, nonsuper- 
radiant sector of the spectrum was shown. The spec- 
trum of massive fields in a Kerr background demonstrates 
a number of phenomena for boundary conditions other 
than quasinormal, such as amplification of an incident 
wave (superradiance) and superradiant instability (see 
Sec.|VIILFJfor details). 

For the D > 6 Tangherlini black holes the fundamental 
mode does not reach quasiresonance, but approaches the 
field mass asymptotically. Therefore, there is another 
interesting effect: the appearance of two concurrent long- 
lived modes with almost the same damping rates. For 
some values of masses of the field /i and of the black hole 
M, the superposition of these two oscillations appears 66 . 

The quasinormal spectrum of a massive vector field in 
the Schwarzschild and Schwarzschild-anti de Sitter back- 
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FIG. 6 Three dominant quasinormal frequencies of spheri- 
cally symmetric massive scalar field perturbations in the back- 
ground of a Schwarzschild black hole (M = 1) for various 
scalar field mass m. 



grounds was studied in 2 ^ 2 -. In the Schwarzschild back- 
ground the fundamental mode of the massive vector field 
behaves qualitatively similar to the massive scalar one: 
As the field mass grows, QNMs reach quasiresonance and 
disappear. The higher overtones become purely imagi- 
nary at sufficiently large values of the field mass. The 
asymptotically high overtones also do not depend on the 
field mass. 

In the asymptotically de Sitter and anti-de Sitter back- 
grounds, quasiresonances do not exist because the mass 
term does not change the boundary conditio n 48 ! 61 . 

Massive fields are short ranged, so their quasinormal 
ringing will almost certainly not be observed in near fu- 
ture experiments. Nevertheless, quasinormal modes of 
massive fields are not only of academic interest: There 
are numerous of situations where massless fields gain an 
effective mass. This concerns a self-interacting mass- 
less scalar field 21 ^, when one takes account of interaction 
with a black hole, black holes in models with small ex- 
tra dimensions 81 and black holes immersed in a magnetic 
fieldZii. 



H. Dependence of the fundamental frequency on black 
hole parameters 

1. Black hole mass: It is possible to fix one scal- 
ing symmetry in the wavelike equation, which al- 
lows us to choose the black hole mass M = 1. In 
order to convert the frequencies calculated in geo- 
metrical units into kHz, one should multiply u) by 
2n(5.U2kHz)M@/M. For example, for the funda- 
mental mode of the Schwarzschild black hole per- 
turbation 1 = 2, 

loM w 0.3737- 0.0890i. 

Then, one can find that for a black hole of 10 solar 
masses the oscillation frequency is 

v k. (1.2074 -0.2875i)fci3*. 

The same dependence on the black hole mass (lj ~ 
M _1 ) was shown for evaporating black holes 21 ^, 
where the dynamics of the dominant quasinormal 
frequency was studied when the black hole mass 
changes. As the black hole mass grows, the real 
and imaginary parts of the frequencies decrease. 
The same was observed for the Vaidya background 
in l 217 ' 218 and /J-dimensions 2 ^. For the evaporat- 
ing three-dimensional AdS black hole the real and 
imaginary parts of the quasinormal modes decrease 
during the black hole evaporation 2 ^. 

2. Electric charge: For small values of the elec- 
tric charge Re(w) and the damping rate grow as 
the charge increases. At some point close to 
the extremal charge Re(u;) reaches its maximum 
value and then decreases. It has been observed 
that this point approximately coincides with the 
point of phase transition, when the heat capac- 
ity of the Reissner-Nordstrom black hole becomes 
positive^ 21 .. The damping rate also reaches its max- 
imum after the critical point. However, most prob- 
ably this is just a numerical coincidence because it 
takes place only for a particular multipole number. 

When we consider charged fields in the background 
of a charged black hole, we must take into account 
electromagnetic interactions between the charge of 
the field and the charge of the black hole. In most 
cases, as the product of the charges grows, the 
real part of the quasinormal frequency grows, while 
the damping rate decreases. This does not hap- 
pen if the black hole charge is close to its extremal 
value. At the near extremal charge the behavior 
of the quasinormal frequencies becomes more com- 
plicated: ui as a function of the black hole charge 
describes spiral-like curves in the complex plane 2 ^. 

3. Rotation parameter: If the angular momentum 
of the particular component of perturbation (m > 
see Sec. IIII.Lf) coincides with the direction of the 
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black hole rotation, the energy is extracted from the 
black hole. This causes decreasing of the damping 
rate and, for some black holes, leads to an insta- 
bility (see Sec. IVIII.C|) . When the direction of the 
angular momentum of a perturbation component is 
opposite to the rotation of the black hole, the per- 
turbation looses additional energy when interacting 
with the black hole. That is why they usually de- 
cay quicker than components with m > 0. The 
actual oscillation frequency grows with m, which 
can be seen, e. g., for the rotating BTZ black hole 
(see Eq. (|4.12|Q . whose QNMs can be obtained as 
an exact solution 202 . For the Kerr black holes we 
observe basically the same behavior—. 

4. Cosmological constant: A positive cosmological 
constant suppresses both the oscillation frequency 
and the damping rate^l for all black holes and all 
types of perturbations 68 ' 103 ' 108 . As the black hole 
approaches its extremal size, i. e., the event horizon 
approaches the cosmological horizon, the quasinor- 
mal frequencies tend to zero, being proportional 
to the surface gravity at the event horizon (see 
Eq. (glUl) ). 

5. Extra dimensions: The radius of the Tangherlini 
black hole horizon r+ is related with its mass M as 



-£>-3 



(D - 2)MA D -2 



4tt 



where Ad-2 is the area of a unit (D — 2)-sphere. 
Therefore, if one compares scattering properties 
around black holes living in worlds of different num- 
ber of space-time dimensions D, it is convenient to 
measure all the quantities in units of the black hole 
radius. 

For a black hole in the scenario with large extra di- 
mensions, the increasing of D leads to increasing of 
both the real frequency and the damping rate. The 
damping rate grows quicker than the oscillation 
frequency and thus, the higher-dimensional black 
hole is a worse oscillator than the lower-dimensional 
one. For the brane-localized fields, the real part of 
the frequency decreases, while the damping rate in- 
creases. 

6. External magnetic field: For the Ernst black 
hole, which is a black hole immersed in an asymp- 
totically uniform magnetic field, it was shown that 
the first order correction from the external mag- 
netic field upon the motion of the scalar field has 
the form of an effective mass A^// ~ m 2 B 2 , where 
m is the azimuthal number and B is the strength of 
the magnetic field^Ll. For a massive charged scalar 
field in the Kerr- Newman background there are two 
main effects on the quasinormal spectrum: the Zee- 
man shift of the particle energy in the magnetic 
field and the difference of values of an electromag- 
netic potential between the horizon and infinity, i.e. 
the Faraday inductior 21 ' 



7. The Gauss-Bonnet term diminishes the damping 
rate of the fundamental mode, implying a black 
hole of the same radius is a better oscillato r 140 ' 144 . 



V. QUANTIZATION OF THE BLACK HOLE AREA 

Even though we do not have a fully noncontradictory 
theory of quantum gravity at hand, there are numerous 
attempts to study black holes using the general princi- 
ples of quantum theory. There are indications that the 
quasinormal spectrum of a black hole possibly allows us 
to find the law of a black hole area quantization. Accord- 
ing to Bohr's correspondence principle, in the regime of 
large quantum numbers, i. e., when the quantum effects 
are small, one must reproduce classical conservative (adi- 
abatic) quantities. For stationary black holes, one such 
quantity is its mass, which must have a discrete spectrum 
in the quantum theory. 

It was suggested in22i that the highly damped oscilla- 
tion frequencies of a black hole are associated with the 
scmiclassical limit of the black hole mass transition, while 
their damping rates correspond to the relaxation time. 
The relaxation time for asymptotically high overtones 
goes to zero, which is compatible with the semiclassical 
approximation. 

Thus, one can write the spacing for the Schwarzschild 
black hole mass as 



AM = HAuj, 



(5.1) 



and, consequently, for the black hole area as 

AA = A(16ttM 2 ) ~ 32ttMAM = 32irhMAu). (5.2) 

S. Hod suggested 2 ^ taking the asymptotical limit of 
the real part of the quasinormal frequency as Aw: 



Alu = He(u!oo) 



In 3 
8ttM ' 



However, some difficulties were found from this conjec- 
ture: 

1. The value of the Re(woo) is non- vanishing for grav- 
itational perturbations, while, e. g., for the vector 
and Dirac perturbations it is zero (see Sec. IIV.Fl) . 

2. The limits of zero rotation and zero charge do not 
commute with the large overtone limi t 226 ' 227 . For 
an infinitesimal rotation the limit of the real part 
and, consequently, the area quantum becomes ar- 
bitrarily small. 

3. Such a value of Auj corresponds to the transitions 
from the ground state of the black hole to a state 
with large "^ 28 - 



All these difficulties can be removed if we consider a 
transition n — > (n — 1) 1 and associate the energy 
with the eigenfrequency of some damped oscillator. This 
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eigenfrequency does not coincide with the real part of the 
quasinormal frequency 2 ^. 

Indeed, consider the damped harmonic oscillator which 
can be described by the wave equation 

£ + 2 7 £ + w& = 0. 

The resonant frequencies of the oscillator are the roots 
of the characteristic equation 

U! 2 + 2lJU! — OJq = 0. 

The solutions to this equation are 



cj = ±yw§ - 7 2 - i 7 . 

Thus, oj 2 = Re(w) 2 + Im(cj) 2 . 

The n — > (n — 1) transition corresponds to the fre- 
quency spacing 

1 



Aw = AVRe(w) 2 +Im(w) 2 ~ Alm(w) 
Then, the area quantum is 



4M 



AA = 8irh = 8nl 2 , 



(5.3) 



where l p is the Plank length. Although the result ob- 
tained was initially demonstrated for the Schwarzschild 
black hole, it was found in a similar way that for the 
Reissner-Nordstrom and Kerr solutions the area quan- 
tum is the sam o 229 ' 230 . Thus, the expression for the area 
quantum (|5.3j) seems to not depend much on the details of 
the black hole geometry. The quantization of the black 
hole area in three and five dimensions was considered 

,231 



Within loop quantum gravity the basis of the Hilbert 
space is given by spin networks. The surface acquires the 
area due to spin edges puncturing 2 ^ 



A = 87rZ 2 7^v/^Tl), 



(5.4) 



where 7 is the Barbero-Immirzi paramete r 233 i 234 , and ji 
is a (half)integer spin. 

It was suggested, that the asymptotic behavior of the 
quasinormal modes allows us to fix 7 by taking^ 

A A — A 

where the minimal surface is given by 

Amin — &7tlp^{ \/ jmin (jmin 1)- 

Here j min is the lowest possible spin. Then, from (|5.3p . 
we find 



T \/ jmin (jmin 1) 1- 

However, in order to reproduce the Bekenstein entropy, 
we need an additional assumption: Either that j = 



1/2 edges puncturing the horizon is suppressed due 
to some unknown dynamics 2 ^ or an introduction of 
supersymmetry 2 ^!. 

As was pointed out in 2 ^, the semiclassical limit re- 
quires consideration of large j, for which the area quan- 
tum is 

AA = 8vr/ 2 7 A V7IJTT) ~ 4^ 2 7 , j > 1. 

This fixes 7 = 2. 

The black hole entropy is given by22&, 



S = ^ = 2,— + Otll 



A 

AA 



(5.5) 



Finally, we note that all the above arguments are 
rather more speculative than mathematically strict. The 
correlation among the quasinormal modes, area quanti- 
zation, and the Barbero-Immirzi parameter, if it exists, 
is probably much more subtle than one could naively ex- 
pect. 



VI. OBSERVATION OF GRAVITATIONAL WAVES 
A. Gravitational waves from black holes 

Black holes can be perturbed due to different processes, 
such as formation of a black after collapse of supernovae, 
black holes, and/or neutron stars mergers, or an infall 
of matter into a supermassive black hole 2 ^ 9 - - - 2 ^ 2 -. One of 
the most frequent phenomenon is the perturbation pro- 
duced by accretion of matter. Unfortunately, the esti- 
mations of the gravitational-wave amplitudes show that 
these waves are too weak to be detected by gravitational- 
wave antennas within the present sensitivity. However, 
within the hyperaccretion scenario gravitational waves 
are potentially detectable for highly spinning black holes 
(a ~ 0.98) by the next-generation detectors LIGO H 2 ^ 3 -, 
LISA 2 ^^, and ET— . 

The estimations for gravitational waves from black 
holes produced by stellar collapse show that it may be 
possible to detect gravitational waves from the collapse 
of stars whose masses are at least 300M(^.. These 
heavy primordial stars have not been observed; however, 
currently there is existing evidence of such intermediate- 
mass black holes, and we can expect that their ringing 
is observed by the next generation of gravitational-wave 
antennas 2 ^. 

In fact, the only realistic scenario for the detection of 
gravitational waves from astrophysical black holes with 
the help of the current antennas is observations of colli- 
sions of black holes and/or neutron stars. Unfortunately, 
the estimated binary merger event rate is also small for 
the present sensitivity (see Sec. IVI.Bj) . Nevertheless, the 
recent estimations show that the black hole-black hole 
inspirals are more likely sources than systems of neutron 
stars and can be detected in the optimistic scenario even 
with the help of the current detectors such as LIGO or 
VIRG0 2 42. 
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The ringing of supermassive black holes in the centers 
of galaxies 2 ^ - can probably be detected by the future 
space-based detector LISA. There are many preliminary 
estimates for event rates of supermassive black hole merg- 
ers, in which the number depends on the model and varies 
from several to thousands of events per yea r 251 ' 252 . LISA 
should also detect gravitational waves from intermediate- 
mass black hole mergers, but the estimates of such events 
are even more uncertain: There are predictions that we 
will be able to observe several events per year^^. 

Since there are a number of works and reviews that dis- 
cuss the possibility of observations of gravitational waves 
from black holes, we do not discuss this topic in more de- 
tail here, but instead refer the reader to a recent review^. 



B. Ringing of neutron stars 

Black holes are not the only source of detectable grav- 
itational waves. Another source of gravitational waves 
are neutron stars, which have been extensively studied 
during the past four decades. Efforts in this direction 
are motivated by advances in the first generation of large- 
scale interferometric gravitational- wave detectors (LIGO, 
GEO600 and Virgo) , which have reached the original de- 
sign sensitivity 2 ^. During the next five years, the second 
generation of gravitational- wave detectors is expected to 
reach their level of sensitivity, which is one order better 
than the first generation. Gravitational-wave observa- 
tions are expected to be complementary to the informa- 
tion obtained from electromagnetic observations 2 ^ 5 . This 
could probably constrain various theoretical models of 
neutron star a 256 ' 257 . 

The most promising sources of gravitational waves 
among various types of stars are as follows: 

• Binary evolution: Binary systems, i. e., systems 
of two stars orbiting around their common center 
of mass, radiate due to rotation and, at the last 
stage, the orbits shrink causing an increase in the 
amplitude of the gravitational radiation. The am- 
plitude of the signal can be calibrated by the two 
masses. From the observed signal one can extract 
the spin rates of the objects and the distance to the 
source 2 ^. The horizon distance for the current an- 
tennas is about SOMpc, which gives the estimated 
frequency of the binary merging of one event per 
25-400 years. Since the second generation detec- 
tors are expected to provide a horizon distance of 
about 300Afpc, which corresponds to several events 
per year, one can expect that gravitational waves 
from binary merging will be detected by the ad- 
vanced LIGO and Virgo detectors 2 ^. 

• Collapse of a supernovae core is described by 
complicated dynamics, which requires solving Ein- 
stein's equations together with equations of state 
for the matter and taking account of neutrino in- 
teractions. Recent numerical simulations of core 



collapse and the evolution of the collapsed remnant 
indicate that gravitational waves can be detected 
by modern gravitational antenna s 260 ' 261 . 

• Rotating stars with nonaxially symmetric defor- 
mation of the star's crust, core, or in an exter- 
nal magnetic field are expected to lose their ro- 
tational energy, which is emitted as gravitational 



waves- 



261,262 



Radio pulsars are promising candi- 
dates for the observation of gravitational waves be- 
cause electromagnetic observations can be used to 
target searches of the gravitational signals. Ac- 
cording to a recent simulation, the deformation of 
the star's crust can be large enough to produce the 
gravitational-wave signal which can be detected by 
advanced LIGO and Virgo detectors 2 ^. 

• Oscillations of stars: Neutron stars have differ- 
ent oscillating modes, which can produce gravita- 
tional waves. These modes could allow us to probe 
internal processes of the stars. Numerical simula- 
tions, based on realistic equations of state, provide 
a set of empirical relations between the frequen- 
cies of the oscillations and star parameters. Some 
of these modes can be detected by the new gener- 
ation of the gravitational-wave detectors. If such 
gravitational waves are detected, we will be able to 
deduce a star's mass and radius with relative error 
of less than l%2SS.2Si. 

Because of their high density, neutron stars have to be 
excellent sources of gravitational waves. The detection of 
their ringing, in particular with the help of the next gen- 
eration of gravitational-wave detectors such as the Ein- 
stein Telescope, will give us the possibility to probe the 
physics of their interior (see topical review 2 ^). 



VII. LATE-TIME TAILS 

A. Generic approach 

As mentioned, as a manifestation of the noncomplete- 
ness of the set of quasinormal modes, at sufficiently late 
time the quasinormal modes are suppressed by the ex- 
ponential or power-law tails. The generic approach to 
analysis of the late-time tails, which we relate here, was 
suggested in 2 ^ 6 -. We start from the wave equation (|2.12|) 
without implying the stationary ansatz ^ ~ e %UJt , i. e., 



d 2 <f d 2 ^ XT/ \ t n 

We define the pair of the following operators: 

d 2 d 2 



D 

D(u>) 



drl 

<£_ 

dr 2 



- uj 2 + V(r.). 



(7.1) 

(7.2) 
(7.3) 
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The wave function ^(V*, t) can be written for t > as 
the integral over the spatial coordinate 

*(r„,t) = J driG(r*,ri;t)d t y(rl,0) + 

J dria t G(r.,ri;t)*(r'.,0), (7.4) 

where G(r*,7^; t) is the retarded Green's function 

IX?(r*,ri;*) = (7.5) 

The Fourier transform reads 

5(w)G(r*,ri;w)=«(r.-ri). (7.6) 

The boundary conditions for black holes and stars differ, 
though the late-time behavior does not because it de- 
pends only upon the asymptotics of the effective potential 
at spatial infinity. Therefore, we think of the radial tor- 
toise coordinate r* as the one that changes from zero to 
infinity < r* < oo (half-line problem) for some nonsin- 
gular systems, such as relativistic stars, or as the one that 
changes from the event horizon to infinity — oo < r* < oo 
(full-line problem) for black holes. The boundary condi- 
tion at infinity for both cases is the same: pure outgoing 
waves at infinity, while at the left boundary the boundary 
conditions read 

9(r*,u) oc e~ iur ', rw-oo (7.7) 

for black holes, and 

$(r*,a;) 0, r* -> (7.8) 

for some nonsingular objects. For simplicity we consider 
potentials which vanish at both boundaries. 

We define the two functions f(ui, r*) and g(w,r*) which 
are solutions to the equation 

D{u)f(u,r.)=D(u)g(u,r.)=0, 

where / satisfies the left boundary condition (at the hori- 
zon or in the origin) and g satisfies the right boundary 
condition (at spatial infinity). The normalization condi- 
tion reads 

lim [e-*" r -0(w,r,)] = l; 
and for the half-line problem 

/(w,r„=0) = 0, />,r,=0) = l, 
while for the full-line problem 

lim [e iur 'f(u,u)] = l. 

r * — v — oo 

The Wronskian has the form 

W(w) = W(g,f)=gd r J-fd r .g 
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FIG. 7 Singularity structure of G(r*, r»; lu) in the lower half ui 
plane, and the contributions of the corresponding residues to 
the Green's function G(r*, r' t \t). The figure is taken from2£&. 



and is independent of r„. Then, one can find that 

G(r.,rl;w) = I 

[f(u,ti)g(w>r*)/W(u) , K < r. 

(7.9) 

Now one needs to integrate / from the left and g from 
the right, until some common point. 
We consider the inverse transformation 

oo 

G(n,r» = J dtG{r*y^ty"\ (7.10) 
o 

and for t > we close the contour in the complex w-plane 
by a semicircle of radius C in the lower halfplane, taking 
the limit C — >• oo. Then, one must take into consideration 
poles of G for Im(w) < (Fig. [7])- Further, one can 
distinguish three different contributions to G. 

1. For a wide class of potentials, for instance, if 
F(r* -> oo) ~ r- a (lnr*)P, [3 = 0, 1, the func- 
tion g(uj, r*) has singularities on the —Imw axis 
and has, thereby, the form of a branch cut along 
the imaginary axis. This branch cut contribution 
to the Green's function is called the tail contribu- 
tion. 

2. At some complex to = Uj, the Wronskian W(oj) has 
zeros and, therefore, the functions / and g are lin- 
early dependent, i. e., one can find a solution that 
satisfies both left and right boundary conditions, 
i. e., uij are quasinormal modes. This contribution 
is naturally called the quasinormal modes contri- 
bution. This contribution decays exponentially as t 
increases, so that for the asymptotically late-time 
behavior t — > oo, the quasinormal modes contribu- 
tion is negligible for most cases (although not for 
AdS space-times, see Sec. IIII.A"|) . 
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3. The contribution from the semicircle \oj\ = C, 
C — > oo corresponds to large lu or short time. This 
contribution vanishes after some time and therefore 
is called the prompt contribution. 

Now we are in a position to summarize the results ob- 
tained for asymptotically late-time behavior of various 
black holes. 



B. Late-time tails of massive and massless fields 

The first study of late-time tails was performed by R. 
Priced! who analyzed the asymptotic behavior of the 
linearized wave equation for the Schwarzschild black hole 
and showed that perturbations of the massless scalar and 
gravitational fields decay according to the power law 



(7.11) 



where i is the multipole number. Soon after, Bicak 
found 2 ^ that the Reissner-Nordstrom tails for the scalar 
field are different for nonextremal and extremal charge, 



1*1 



t -(2l+2)^ iiQ<M . 

trV+ 2 \ if Q = M. 



(7.12) 



The asymptotic tails for the Kerr space-time were shown 
to be the same as for the Schwarzschild one. We note 
that the above results hold when the problem allows for 
"good" initial conditions of compact support. 

Unlike asymptotically flat solutions, de Sitter 
black holes have qualitatively different late-time 
behavio r 269 ' 270 . First, the asymptotic tails are not power 
law but exponential, and second, the I = mode does 
not go to zero, as takes place for the Schwarzschild case, 
but asymptote to some constant 



1*1 



-lk a t 



if 1=1,2, 



1*1 



|*o| + |*i|e 



-2k c t 



if e = o, 



(7.13) 
(7.14) 



where k c is the surface gravity on the cosmological hori- 
zon. 

The late-time behavior of massive fields is drastically 
different from massless ones. The late-time tails are not 
purely decaying but intensively oscillating for massive 
fields (see, for instance, Fig. UTI) . One of the reasons is 
that for massive fields tails appear already for Minkowski 
space-time, which serves as a dispersive medium for mas- 
sive scalar field. Using the exact Green's functions it was 
showr>2Zi that the late-time tails of the scalar field in flat 
space-time are 



|*/iat| ~ t 1 5 sin(^t), 



(7.15) 



where fi is the inverse Compton wavelength. The 
Minkowski space-time tail shows itself in the black hole 
tails at the so-called intermediately late time, i. e., when 

1 < t/M < (fiM)- 3 , 



for massive scalar and Dirac fields. For the Proca field, 
the intermediate late-time tail depends on the polar- 
ization of the field and can be either i _ ( f+1 / 2 ) sin(jUt), 

t -(i+3/2) sm ( Mt)) Qr t-it+6/2) sin (^)272. 

In the presence of a black hole the massive scalar field 
decays at asymptotically late time (which now can be 
defined as t/M > (fiM)~ 3 ) as 2 -^ 



|*| ~t e sin(/zt), t — > oo, 



(7.16) 



independently of the angular number I. This asymptotic 
decay is quite universal not only because of independence 
on £ but also because it takes place for a quite gen- 
eral class of black holes, which includes Schwarzschild^ 7 - 3 -, 
K err- 27i or dilator* 2 ^ solutions as well as for other massive 
fields (Dirao22£ and Proca 2 ^ 2 -). One may expect, there- 
fore, that this behavior (|T. 16[) is universal for all massive 
fields and does not depend much on the details of the 
black hole horizon geometry^ 2 -. Although this supposi- 
tion sounds too strong and was not checked until recently. 
Finally, it should be noted that the Price decay law 

|ljr| ~ £-(2^+3) ho l ds for 

some time-independent back- 
grounds, while for timc-dcpcndcnt solutions a different 
law |*| ~ £~( 2£ + 2 ) was recently found 2 ^ 6 -. 



C. Are asymptotic tails pure nonlinear phenomena? 

In the linear approximation, the study of asymptotic 
tails in higher than four dimensions was initiated by Car- 
doso et al. ir*22I. There it was stated that the late-time 
decay law reads 



I* I ~ f -(2e+3D-S) 
I 1 even \ L 



(7.17) 



for an even number of space-time dimensions D > 4 and 
is 



(7.18) 



for odd D=is. Now the case D = 4 looks quite special. 
Further, it was shown that in odd dimensions the late- 
time tails are independent on the details of the black hole 
geometry and exist already in the Minkowski space-time 
even for massless fields. Therefore, it did not come as a 
surprise that odd-dimensional Gauss-Bonnet black holes 
have the same decay law at asymptotically late timei 3 -^. 

Recently, however, the above results, obtained in the 
linear approximation for D > 4 black holes and proba- 
bly for some D = 4 solutions, have been challenged inSZ 9 . 
by consideration of the fully nonlinear perturbations of 
the Z)-dimensional Schwarzschild black hole. InSZ 9 . it is 
stated that the nonlinear approach gives slower decay at 
asymptotically late time for all D, except D = 4, so that 
the asymptotic tails seem to be essentially a nonlinear 
phenomenon. This implies that the D = 4 case appears 
as an exception when both linear and nonlinear analyses 
lead to the same result. The possible nonlinear effect on 
asymptotic tails was also reported ir*22°-, although, non- 
linear analysis of black hole tails is still poor and needs 
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further effort in order to prove or disclaim the results 
obtained in the linear approximation. 



VIII. STABILITY 

A. Gravitational (in)stability: General aspects 

Gravitational stability of any static or stationary so- 
lution describing a compact object which is expected to 
exist in nature is a principal issue of the physical ade- 
quateness of the solution. Objects which are unstable 
under small perturbations will be inevitably destroyed 
by them and, thereby, simply cannot exist. A compact 
object can be a black hole or brane, a star or a worm- 
hole, or even more exotic objects such as a black ring, 
black saturn, or a naked singularity. Any of them can be 
accepted as properly described by its metric only after 
the metrics is proved to be stable. 

Since dynamical perturbations of a gravitational sys- 
tem can usually be described by a single wavelike equa- 
tion (12.12)) (or a set of such equations) , one can analyze 
stability against linearized perturbations in the follow- 
ing way. Under existence of the well-posed initial value 
problem, the complete spectrum of quasinormal modes, 
we judge about stability: if all quasinormal modes are de- 
caying in time, the solution is evidently stable. If there 
is at least one growing mode, the space-time is unstable 
with the instability growth rate, which is proportional to 
the imaginary part of the growing QNM. 

If the effective potential V e ff in the wave equation 
(I2.12j) is positive definite, the differential operator 



is a positive self-adjoint operator in the Hilbert space of 
square integrable functions L 2 (r*, e?r»). Then, there are 
no negative (growing) mode solutions that are normaliz- 
able, i. e., for a well-behaved initial data (smooth data 
of compact support), all solutions are bounded all of the 
time. 

If the effective potential is not positive definite ev- 
erywhere outside the horizon, a special trick, called the 
Friedrichs extension, can be used, which allows us to 
prove stability (if there is one) in some cases. If the 
effective potential V e f / is bounded from below, the oper- 
ator A is a symmetric semibounded operator for smooth 
functions with compact support in Then, A can be 
extended to a semibounded self-adjoint operator in such 
a way that the lower bound of the spectrum of the ex- 
tension is the same as that of A. Such an extension is 
unique for asymptotically flat space-times, because the 
range of r* is complete (—00, +00), and thereby the time 
evolution of the wave function is uniquely determined 
by the chosen initial data. Finally, if one can find the 
positive definite extended operator, this means that the 
space-time is stable. 



For asymptotically anti-de Sitter space-times, the 
range of r*, which is (— 00, r°), where r° is some con- 
stant, is not complete and the Schwarzschild wedge 
is not globally hyperbolic. Nevertheless, upon impos- 
ing the Dirichlet boundary conditions (|3.3|) . one can 
use the Friedrichs extension as well^r— . Unfortu- 
nately in the majority of interesting cases, such as 
higher-dimensional Schwarzschild- AdS, Schwarzschild- 
dS, Reissner-Nordstrom black holes and various other 
generalizations, the Friedrichs extension method (which 
is also called the S'-deformation method) does not work 
well, because it is difficult to find an appropriate ansatz 
for the function which performs the extension of the oper- 
ator. Therefore a practical tool for testing stability in all 
those cases is numerical investigation of the quasinormal 
spectra 68 i 127 . 



B. Gravitational stability in four dimensions 

The first work on the black hole stability (by Regge 
and Wheeler—) proved stability of the (3+l)-dimensional 
Schwarzschild black holes against axial type of perturba- 
tions (though polar perturbations were considered there 
as well but with a mistake which was corrected later by 
Zeriili 2 ^!). Since that time, the Schwarzschild solution 
gained status of a physically adequate model for neutral 
isolated nonrotating black holes. 

The next step was an analysis of perturbations and 
proof of stability of Reissner-Nordstrom black holes per- 
formed, almost at the same time and independently, by 
N. Sibgatulin and G. Alexee-v 2 ^ and by Moncrief 2 ^. 
In222 in addition to the decoupling of variable and deduc- 
tion of the master wave equations, the gravitational scat- 
tering around Reissner-Nordstrom black holes was con- 
sidered that included construction of the Green's function 
for the wave equations. 

The Kerr black hole was shown to be stable against 
linear perturbations by Teukolsky and Press by using the 
first-order WKB approac h 284 ' 285 . A mathematical proof 
of stability was performed much later by Whiting 2 ^. 

The (3 + 1) Schwarzschild-de Sitter and Reissner- 
Nordstrom de Sitter black holes were proven to be stable 
as well2£, and the stability of the Schwarzschild-anti-de 
Sitter solution was shown mi 5 -^. The stability of Kerr- 
de Sitter black holes was shown recently ir>22£, and of 
Kerr- AdS black holes ini^. 

Finally, the D = 4 string theory inspired black holes, 
such as dilaton black holes, Born-Infeld black holes, and 
Gauss-Bonnet coupled to dilaton black holes, were shown 
to be gravitationally stable as wel l 92 ' 288 " - 2 ^ 2 -. Although 
the stability analysis of the Born-Infeld black holes^ 8 - 8 - is 
still incomplete and includes only axial perturbations, 
leaving therefore a chance of instability in the polar 
modes. 

All of the above mentioned works are summarized in 
Table Hill We concluded, therefore, that all of the consid- 
ered here four-dimensional black holes, tested for stabil- 
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TABLE III Stability of four-dimensional black holes. All of 
the four-dimensional black holes considered here are stable, 
except Kerr-Newman black holes and its string theory gen- 
eralization (dilaton-axion black holes) for which the variables 
cannot be decoupled in an easy way. 



Black hole solution (parameters) 


Publication 


Schwarzschild (M) 




Reissner-Nordstrom (M, Q) 


2o2 , zoo 


Schwarzschild-dS (M, A > 0) 


96 


Schwarzschild-AdS (M, A < 0) 




Reissner-Nordstrom-dS (M, Q, A) 


96 


Kerr (M, J) 


2o4, ZOO 


Kerr-dS (M, J, A > 0) 




Kerr-AdS (M, J, A < 0) 


163 


Kerr-Newmann (M, J, Q) 


? 


Kerr-Newman-A(dS) (M, J, Q, A) 


? 


Dilaton (M, Q, <f)) 


92.2S9.292 


Dilaton-axion (M, Q, J, (j>, ip) 


? 


Dilaton-GB (M, 0) 


290,291 


Born-Infeld (M, Q) 


axial*** 


Black universes (M, (f>) 


293 


BHs in the Chern-Simons theory (M, /?) 


294 



ity, proved to be stable. Although Kerr-Newman black 
holes, and their string theory generalizations, which in- 
clude axion and dilaton fields, are still not tested for sta- 
bility because the perturbation equations do not allow 
for an easy decoupling of the angular variables. 



C. Gravitational instabilities in higher than four dimensions: 
Gregory-Laflamme and non-Gregory-Laflamme instabilities 

We have seen in the previous subsection that four- 
dimensional black holes are usually stable against gravi- 
tational perturbations. The situation in higher dimen- 
sions is much richer, where one has various instabili- 
ties. A wide class of D > 4 objects, black strings and 
their various generalizations, such as black branes, suf- 
fer from a general type of gravitational instability, called 
the Gregory-Laflamme instability 2 ^-. The essence of this 
phenomenon can be easily understood when considering 
linear perturbations of black strings 296 . 

Unlike Kaluza-Klein black holes, the black string met- 
ric is a solution to the Einstein equations in five- or 
higher-dimensional gravity which has a factorized form 
consisting of the Tangherlini black hole and an extra flat 
compact dimension. According to the brane-world sce- 
narios, if the matter localized on the brane undergoes 
gravitational collapse, a black hole with the horizon ex- 
tended to the transverse extra direction will form. This 
object looks similar to a black hole on the brane, but is, 
in fact, a black string in the full D-dimensional theory. 



The generic black string metric has the following form: 

ds 2 = g ab dy a dy b + dz 2 , (8.2) 

where g a b describes the black hole behavior on the brane, 
while the "extra" dimension is in the z-direction. The z- 
direction is periodically identified by the relation z = 
z + 2ttR. 

The perturbation of the Einstein equations for black 
strings can be reduced to scalar, tensor, and vector types 
of perturbations in the same way as was done for higher- 
dimensional black holes (see Sec. HI|) . An analysis of 
H. Kudoh222 showed that vector and tensor gravitational 
perturbations are stable, as well as I = 1, 2, ... scalar per- 
turbations. The only unstable type of perturbations is 
£ = scalar gravitational type, also called an s-wave. 
This s-wave shows the Gregory-Laflamme instability at 
long wavelengths in the z-direction. 

The zero mode (k z = 0) of the scalar gravitational per- 
turbation with I = corresponds to a shift of the mass of 
higher-dimensional Schwarzschild black holes, and, there- 
fore, this mode is not dynamical due to BirkhofFs the- 
orem. However, the Kaluza-Klein mode with k z =/= 
and £ = is essentially different from the gravitational 
perturbations of Schwarzschild black holes, because it 
cannot be interpreted as an infinitesimal change of mass 
of the black string, i. e., it corresponds to a dynamical 
perturbation. Thus, from the viewpoint of an effective 
theory on the z — const plane, the Gregory-Laflamme 
instability is apparently related to the inapplicability of 
BirkhofFs theorem. An interesting membrane illustra- 
tion of the Gregory-Laflamme instability was suggested 
in228. 

A qualitatively different type of instability, which oc- 
curs not at the s-wave but at the first multipoles, takes 
place for some higher-dimensional black holes. For ex- 
ample, the _D-dimensional Reissner-Nordstrom-de Sitter 
black hole, described by the metric 

ds 2 = f(r)dt 2 - f-\r)dr 2 - r 2 da 2 D _ 2 , (8.3) 

where do~ 2 D 2 is the line element of a unit (D — 2)-sphere, 

f n = — — ^ 2Ar 2 

JV> r D-3 + r 2(D-3) (D_ !)(£>_ 2)' 

is unstable for large values of the charge and cosmological 
termi 2 ^. This instability, unlike the Gregory-Laflamme 
one, occurs at the £ = 2 multipole, and thereby is not 
related to the inapplicability of the BirkhofFs theorem. 
The parametric region of instability can be seen in Fig. [FJ 
while the form of the slightly deformed black hole at the 
threshold of instability is shown in Fig. |H1 

In order to distinguish the latter instability from 
the Gregory-Laflamme one, we call this instability 
a non-Gregory-Laflamme instability. In addition to 
black string, black brane, and Reissner-Nordstrom-de 
Sitter instabilities, there is an instability of higher- 
dimensional black holes in the Einstein-Gauss-Bonnet 
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TABLE IV (In)stability of higher-dimensional black holes 



TV 1111 1 1 ■ / J \ 

Black hole solution (parameters) 


Publication 


Schwarzschild (M) 


Stable for all 


Reissner-Nordstrom (M, Q) 


n i ii r t~~\ k y» i i i i ii 127 

stable lor D = 5,6, ... ,11 and nonextremal charge"* 


ni i*ii in /"\ t A\ /ay * ^ /-\\ 

Schwarzschild-dS (M, A), (M, A > 0) 


f~1 i 11 P TV K /"» 11 fiK 

Stable for £) = 5, 6, . . . , 11— 


f~1 1 1 "1 1 A in /A T A \ /A T A ^ 

Schwarzschild-AdS (M, A) (M, A < 0) 


r~i i ii ■ i — > ii a - j i r ta i-* j-i -i 1 3fl2 

Stable m EM theory for D = 5, 6, . . . , lialia 


Ti ■ "\T 1 i " in /A A" A \ /A A" /^V A r\ \ 

Reissner-Nordstrom-dS (M, A) (M, Q, A > 0) 


T T i ll f* TV t—rd 1 h 197 

Unstable for D = 7, 8, . . . , 114^ 


TV ■ TVT "\ i " A If /TV iT A \ /TV X /~\ A ^ r»\ 

Reissner-Nordstrom-AdS (M, A) (M, Q, A < 0) 


c~i i ii ■ tth /Tii 302 i iii' - j6 

Stable m EM theory^ 4 *- and unstable in supergravity 


Gauss-Bonnet (M, a) 




Myers-Perry and its generalizations (M, J) 




Dilaton (M, Q, </>) 


? 


Dilaton-axion (M, Q, J, <j>, ip) 


? 


Dilaton-Gauss-Bonnet (M, <j>, a ) 


7 



q 




FIG. 8 The parametric region of instability in the right upper 
corner of the square in the p — q "coordinates" for D — 7 (top, 
black), D = 8 (blue), D = 9 (green), D = 10 (red), D = 11 
(bottom, magenta). The units r+ — 1 are used; p = r+/r c — 
l/r c < 1, r c is the cosmological horizon. The charge can be 
normalized by its extremal quantity q = Q/Qext < 1. 
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FIG. 9 The equatorial plane of the black hole horizon hyper- 
surfaces at the edge of stability. The dashed line corresponds 
to the unperturbed BH of unit horizon. The blue, green and 
red lines correspond to the perturbed BHs (for different values 
of A) after decay of all the dynamical modes. 



theor y 140 i 299 ~ ^21. The summary of (in)stability analysis 
in higher dimensions is given in Table llVl 

The case of the Reissner-Nordstrom-AdS solution 
(given by the metric (|8.3j) with a negative cosmological 
constant) deserves special comments. Solution (I8.3|) sat- 
isfies not only the ordinary Einstein-Maxwell equations, 
but also the M — 8 supergravity equations which include, 
among a number of fields, the dilaton. The stability prop- 
erties of the above solution naturally depend upon the 
theory in which the solution is considered. In£ it was 
shown that in the TV — 8 supergravity, large Reissner- 
Nordstrom-AdS black holes are unstable for large val- 
ues of charge, while in the Einstein-Maxwell theory these 
black holes are stable 3 -^. If one remembers that in the 
M = 8 supergravity a dilaton allows for a dynamical 
£ = (i. e., pure spherically symmetric) mode, then it is 
evident that the instability is of the Gregory-Laflamme 
type. This type of instability is certainly impossible in 
the pure Einstein-Maxwell theory, where the Reissner- 
Nordstrom-AdS black hole has a nondynamical s-wave 
mode 3 -^. 

Although the stabilities of nonrotating solutions are 
studied relatively well, the stability of rotating black 
holes and branes in higher dimensions is almost a black 
spot. The most appealing problem is certainly stabil- 
ity of higher-dimensional Myers- Perry black holes and of 
their generalizations. There are two qualitatively differ- 
ent cases here: If all of the angular momenta are equal 
ai = ci2 = ■•• = ajv, then, as in the four-dimensional 
case, the momentum of the black hole is bounded. The 
perturbation equations in this case allow for a particu- 
lar separation of variable a 303 i 304 i 306 . Thus, in the D = 5 
case with a\ = a.2 the separation of variables is possible 
for the so-called zero mode, i. e., for lowest eigenvalues 
of the momentum^. Tensor gravitational perturbations 
allow for the separation of variables for any D and all the 
eigenvalues in this case. 

The second, less symmetric case, when at least one of 
the angular momenta is different from the others, is much 
more interesting. Then there is no bound for the rotation 
parameter a (for D > 5) and an infinitely high angular 
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FIG. 10 The qualitative phase diagram for the black objects 
in D > 6, proposed m2£!£. The horizontal and vertical axes 
correspond, respectively, to the spin and area of a black ob- 
ject. If thermal equilibrium is not imposed, multirings are 
possible in the upper region of the diagram. 



speed of rotation is possible for the Myers-Perry solution, 
at least formally. Intuitively we expect that such a highly 
rotating black hole must be unstable when the centrifugal 
forces exceed the gravitational attraction (see Fig. [T0|) . 
However, until recently, no instability was proved for such 
black holes, because of the impossibility of the separa- 
tion of variable in the perturbation equations. In the 
linearized theory some particular results were obtained 
for tensor gravitational perturbation o 131 ' 132 ' 305 , where no 
instability was found even for asymptotically high rota- 
tion. This apparently means that the instability is in 
the scalar gravitational modes, which are responsible for 
deformation of the horizon. Recently, a static limit of 
perturbations was considere d 308 ^ 309 and the onset of the 
bifurcation point toward a new solution was shown for 
some critical value of the angular momentum. Finally, 
the instability was detected by numerical simulations in 
the fully nonlinear theory for large angular momenta* 3 -^. 



D. Two types of developing of gravitational instabilities in 
time domain 

The development of instability in the time domain 
is different for static, spherically symmetric nonrotat- 
ing black holes and rotating black holes. The spheri- 
cally symmetric black holes have damped QNMs in the 
form of damped oscillations, i. e., with nonvanishing 
real and imaginary parts of ui, while growing, unstable 
modes must be pure imaginary, i.e., nonoscillatory. In- 
deed, we multiply Eq. (|2.12[) by the complex conjugated 
function "J* and assume that the dependence on time 
is ^>(t,r*) = e _iwi \P(r*). The integral of the obtained 
equation reads 



/ = 



d 2 V(r+) 
drl 



+ L0 2 \V(r+)\ 2 -V\V{r*)\ 2 dr+ 




t 



500 



FIG. 11 Time-domain profiles of black string perturbations 
for n = 1 k = 0.84 (magenta, top), k = 0.87 (red), k = 0.88 
(orange), k = 0.9 (green), k — 1.1 (blue, bottom). We can see 
two concurrent modes: for large k the oscillating one domi- 
nates , while near the critical value of k the dominant mode 
does not oscillate (looks like exponential tail), and for unsta- 
ble values of k the dominant mode grows. The plot is loga- 
rithmic, so that straight lines correspond to an exponential 
decay. 



Integration of the first term by parts gives 

oo 

,d*(r*) 



I = **(r*)- 



dr+ 



+ / [u 2 \*(r+)\ 2 -V\*(r+)\ 2 - 



dr± 



dr± = 0. 



Taking account of the boundary conditions (|3.52|) . one 
can find the imaginary part of the integral 

Im(J) = Re(cj)|#(oo)| 2 + Re(w)|#(-oo)| 2 + 

OO 

+2Re(w)Im(w) J |#(r*)| 2 dr* = 0. 



The nonzero real part of the quasinormal frequency im- 
plies that the imaginary part is negative. Therefore, the 
unstable modes (Im(cj) > 0) must have zero real part. In 
other words, for static spherically symmetric black holes 
unstable modes cannot be oscillating. 

In addition to the nonoscillatory character of the un- 
stable modes, there is another distinction of the evolu- 
tion of instabilities in the time domain. The first type 
of the time-domain instability develops immediately af- 
ter the initial outburst. This type of instability is shown 
in Fig. [11] for the black string 3 -!! and in Fig. [12] for the 
Reissner-Nordstrom black holea!2I. 

The second type of time-domain evolution is much 
more exotic: The instability develops after a rather long 
period of damped quasinormal oscillations. This type 
of instability takes place for black holes in the Einstein- 
Gauss-Bonnet theory (see Fig fTB]) 140 ' 299 . The instability 



35 




FIG. 12 Time-domain profile of near extremal RN BH per- 
turbation (D = 11, p = 0.8), q=0.4 (brown), q=0.5 (blue), 
q=0.6 (green), q=0.7 (orange), q=0.8 (red), and q=0.9 (ma- 
genta). The smaller q is, the slower the growth of the profile 
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FIG. 14 Instability region for RNAdS black holes in M = 8 
supergravity. The plot was proposed irA 
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FIG. 13 The picture of instability, developing at large mul- 
tipole numbers: D = 6, I = 8 (red, bottom line), £ = 12 
(green), I = 16 (blue, top line), a = 1.4. Tensor type of 
gravitational perturbations. 



occurs at higher multipoles i, while the first few lowest 
multipoles are stable^. We believe that the stability 
of lowest multipoles gives rise to a prolonged period of 
damped oscillations and the "delayed" appearance of the 
instability in this case. 



E. Correlation between gravitational and thermodynamic 
instabilities 



the Gregory-Laflamme type and, thus, directly connected 
to the s-wave mode and inapplicability of Birkhoff 's the- 
orem. The parametric regions of gravitational and ther- 
modynamic stability do not coincide exactly as shown 
in Fig. PHP ; the region of gravitational stability is larger 
than that of the thermodynamic one and includes the lat- 
ter completely. In other words, the Reissner-Nordstrom- 
AdS black hole that is thermodynamically stable is also 
gravitationally stable but not otherwise, the gravitation- 
ally stable black hole is not necessarily thermodynami- 
cally stable. For large highly charged anti-de Sitter black 
holes the regions of gravitational and thermodynamic in- 
stabilities asymptotically coincide (Fig. Il~4"]) . Most prob- 
ably this happens because in this limit the AdS black 
hole approaches the black brane regime. For the black 
brane, Gubser and Mitra claimed that thermodynamical 
stability is the sufficient and necessary condition for the 
absence of gravitational instability. They also claimed 
that, in the limit of large AdS black holes, gravitational 
and thermodynamic stabilities coincide, while the small 
discrepancy observed in Fig. [14] is probably a numerical 
error—. The physical explanation of the instability con- 
sists of an energetically more preferable state of a set of 
black holes other than of a single black brane: The en- 
tropy of an array of black holes is larger than the entropy 
of the uniform black brane of the same total mass. 



According to Gubser and Mitral, gravitational insta- 
bility is correlated to the local thermodynamic insta- 
bility for the RNAdS black holes in Af = 8 supergrav- 
ity. Mathematically the local thermodynamic instability 
means that the Hessian matrix of the second derivatives 
of the mass with respect to the entropy (and the con- 
served charges or angular momenta) has a negative eigen- 
value. As we have already learned, this instability is of 



F. Superradiant instability 

Consider the classical scattering problem for a massless 
scalar field in the background of the Kerr black hole. The 
problem can be reduced to the wavelike equation (|2.12[) 
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with the following effective potential 

2 \2 



V" 



am 



A, 



(r 2 + a 2 ) 2 
d 1 



A, 



(r 2 + a 2 ) 3 / 2 dr \ r dr ^r 2 + a 2 
r 2 + a 2 - 2Mr, 



, (8.4) 



where A is the separation constant 21 . The effective po- 
tential has the following asymptotic behavior: 



V — > oj 2 , r — > co, 
V — > (oj — milh) 2 , r — > r + . 



(8.5) 
(8.6) 



The asymptotic forms of the solutions near the black hole 
and at spatial infinity are 



*(r„) 



*(r») = T e 



r — > oo, 
7' — > r_i-. 



(8.7) 
(8.8) 



Here 1Z is called the amplitude of the reflected wave or 
the reflection coefficient, and T is the transmission coeffi- 
cient. Since the effective potential is real, the Wronskian 
of the complex conjugate solutions ^(r*) and ^*(r*) 
obeys the relation 

%4-W =i4~ f*?--**?^) =0- (8.9) 

dr* dr* \ dr* dr* J 

Integrating the above relation from r+ until spatial in- 
finity and using the asymptotic form (|8.7p at infinity, we 
obtain 

\K\ 2 = 1 + {i/2u))W\ r+ . (8.10) 
Using the asymptotic form (|8.8p . we find that 



\KV = 1 



- 1 \n 



If \K\ > 1, that is, 



> 1, 

U) 



ill) 



.12) 



the reflected wave has larger amplitude than the incident 
one. This amplification of the incident wave is called the 
superradiance and was first predicted by Zeldovich 3 ^. 
The superradiance effect for Kerr black holes was first 
calculated by Starobinsk y 313 i 314 , who showed that the 
superradiance is much stronger for the gravitational field 
than for scalar and electromagnetic fields. The process 
of superradiant amplification occurs due to extraction 
of rotational energy from a black hole, and, therefore, 
it happens only for modes with positive values of az- 
imuthal number m that corresponds to "corotation" with 
a black hole. The superradiance is absent for fermion 

fieldg 315 1 316 

The condition (|8.12l) can be generalized to perturba- 
tions of fields of other integer spin, as well as to other 



than rotating black hole systems^IZ, for instance, to an 
electrically conductive rotating cylinder or another sys- 
tem with absorption. In the general case the effective 
potential has a nonvanishing imaginary part 

V(r*) = U(r*)+iT(r*). (8.13) 

Then, Eq. (gU) reads 

ijLw(V,y*) = 2r\y(r*)\ 2 , 

and Eq. (|8.10j) gets an extra term 



\TZ\ 2 = 1 + (i/2u)W\ r+ - - / r(r.)*(r.)d(r,). (8.14) 



The condition of superradiance is 



(i/2u)W\, 



1 



T(r*)^(r*)d(r*) > 0. 



The existence of an ergoregion is not sufficient for su- 
perradiance: An appropriate boundary condition is also 
necessary, which, for rotating black holes, is the require- 
ment of only the ingoing group velocity waves at the event 
horizon 3 ^. 

When one has the massive scalar or other integer spin 
field in asymptotically flat or de Sitter space-time or 
massless field but an AdS boundary at spatial infinity, 
the superradiance has unstable modes. The essence of 
this phenomenon can be understood from the plot of the 
effective potentials for massive scalar field in asymptoti- 
cally flat black holes and for massless fields in AdS black 
holes (Fig. I15[) . The effective potentials for both cases, in 
addition to the local maximum, have an extra local min- 
imum far from the black hole which creates a secondary 
reflection of the wave reflected from the potential bar- 
rier. This secondary reflected wave, when incident on a 
potential barrier, will be reflected again at the far region. 
As each reflection from the potential barrier in the super- 
radiant regime increases the amplitude of the wave, the 
process of reflections will continue with increased energies 
of waves and, thus, one has an instability. The superra- 
diance for various black holes was considered in^iS - ^ 2 ^. 

We now discuss in more detail the instability due to 
the massive term. This case implies that at least in four 
space-time dimensions the meta-stable bound states can 
be formed 3 ^ in the valley of the local minimum. These 
bound states are characterized by the product of the field 
mass n and the black hole mass M, i. e., by the ratio of 
the characteristic size of a black hole to the Compton 
wavelength of the particle 



fjtM = GMfi/hc- r + /X c 



(8.15) 



Detweilep222 first estimated the instability growth rate 
for the I = m = 1 bound state in the limit fj,M -C 1 , 



r w 24(a/M)~ 1 (LiM)- 9 (GM/c 3 ), 



(8.16) 
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FIG. 15 The qualitative behavior of the effective potential for massive fields in the Kerr background (left) and for the massless 
fields in the Kerr- AdS background (right). 



where r is an e-folding time. Then, Zouros and 
Eardley 328 computed the WKB instability growth rate 
in the opposite limit (iM ^> 1, 

r a lOV- 84 ^ (GM/c 3 ). (8.17) 

Finally Dolan 3 ^ analyzed the instability for the whole 
range of values of /iM and found that the maximal in- 
stability is for the t = 1, m = 1 state, for /iM < 0.42, 
a = 0.99 and equals t" 1 rs 1.5 • lO' 7 (GM/c 3 )' 1 . it 
was shown that the presence of the strong magnetic field 
can enhance the superradiant instability. The superradi- 
ant instability of massive particles in the vicinity of black 
holes appears always as a negligible process: The reason 
is that the instability growth rate is always small for the 
standard model particles in comparison to the decay rate 
of particles or of the Hawking evaporation rate of black 
holes. 

The superradiant instability due to the presence of the 
AdS boundary was studied i n 331 i 332 for the scalar field 
and i n 132 i 304 for gravitational perturbations. Ini^ 2 - it was 
shown that the instability of the tensor gravitational per- 
turbations of simply rotating D > 6 Myers-Perry black 
holes is due to superradiant modes only and also has a 
small growth rate of order 10~ 12 (r+) -1 . As the instabil- 
ity occurs for small AdS black holes (one can see the plot 
of the parametric range of instability in Fig. I16[) . it is 
certainly negligible and will be suppressed by the violent 
Hawking evaporation. 

Note that at the quantum level superradiant instability 
of a massive particle simply means that the particle is 
made to leave a given superradiant state in favor of a 
"normal" (nonsuperradiant) one. 

IX. ADS/CFT INTERPRETATION OF QNMS 

A. A brief overview of AdS/CFT 

While string theory is often regarded as a theory of 
quantum gravity and grand unification, it also encom- 
passes a wide range of dualities, including AdS/CFT 



pure adS limit 




small bh^ ► large bh 



FIG. 16 The stable region in the parameter plane for the 
simply rotating higher-dimensional asymptotically AdS black 
hole. The figure is taken fromi^. 



which relates strongly coupled field theories to weakly 
coupled gravitational dualsA. As a strong/weak coupling 
duality, AdS / CFT has been applied to the study of the 
strongly interacting quark-gluon plasma as well as to con- 
densed matter systems such as high T c superconductors. 
Here we give some basic ideas about AdS/CFT, and refer 
the reader to a review^ for further study. 

The key theoretical problem of modern particle physics 
is the search for an adequate description of interactions 
of quarks and gluons, i. e., quantum chromodynamics 
that would be valid at all energies. For small values of 
the coupling constant g, the Yang-Mills theory with an 
SU(3) group of symmetry is a good approximation, where 
3 stands for the three colors of quarks. However, at small 
energies and larger distances the coupling constant g is 
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growing, which does not allow us to use an expansion in 
g. An approach that could hopefully solve this problem 
is based on the consideration of the SU(iV) Yang-Mills 
theory in the limit N — > oo^^ 3 ., a theory with an "infinite 
number of quarks' colors" . In this limit the perturbation 
theory becomes much simpler so that all nonvanishing 
diagrams look similar to triangulations of a spher o . It 
is then natural to expect that in the limit N — > oo, the 
Yang-Mills theory can be described within the formalism 
of two-dimensional strings 3 - 3 ^. The contributions to these 
triangulations are expansions in terms of g 2 N , so that 
one keeps g 2 N finite in the regime of large N. 

The so-called conformal-invariant supersymmetric 
Yang-Mills theories have gained a great success in this 
way during the past decade^. According to string theory, 
the world surface of a string has an infinite tower of quan- 
tum excitations which look similar to particles. A finite 
number of these quantum states of a string corresponds 
to massless particles, and the infinite "rest" of the states 
represents massive particles. The mass of a particle is 
proportional to the string tension, so that at large dis- 
tances (r » I s ) only massless particles survive. Among 
massless excitations of closed strings there is one corre- 
sponding to the spin-2 field, that is, to a graviton. Con- 
sequently from string theory at large distances, one can 
deduce the Einstcin-Hilbert theory of gravity, and from 
superstrings - supergravity. The full theory which in- 
cludes gravitation and interacting "matter" fields can be 
described by the system of both closed and open strings 
and a set of the so-called Z?-branes. The Z?-branes are 
submanifolds to whose ends open strings are attached, 
while closed strings can propagate in the bulk as well. 
Our world is supposed to be such a D3-brane (here 3 
is for three spatial dimensions) to which strings are at- 
tached and fluctuate in the bulk. Z?-branes must keep 
some part of the supersymmetry and this stipulates the 
particular form of the space-time geometry in the vicin- 
ity of D-branes. The geometry near the brane's throat is 
the anti-de Sitter one and the anti-de Sitter radius R of 
the .D3 brane is 

R = Ug 2 Nf/\ (9.1) 

where l s is the string length. If a closed string has energy 
which is less than the brane curvature, the string cannot 
overwhelm the brane's gravitational attraction and will 
leave the brane, staying near its throat. If g 2 N > 1, 
one has the regime of a classical superstring in the back- 
ground of the Z?3-brane. 

B. The AdS/CFT vocabulary and interpretation of 
quasinormal modes 

Now we are in a position to relate, briefly, the famous 
AdS/CFT correspondence, the duality between string 
theory and field theory^. On the string theory side we 
have the following constants: the anti-de Sitter radius R, 
string's length £ s and the string coupling g s . On the field 



theory side one has the 't Hooft coupling Ng 2 . The du- 
ality between string and field theory works through the 
following relations between the coupling constants: 

4ng s =g 2 , g 2 N=^. (9.2) 

The main advantage of this mapping is that one can de- 
scribe the strong coupling regime in field theory (large 
g 2 N) by string theory in the regime R 3> is, that is, by 
supergravity. The AdS/CFT vocabulary says that there 
is a correspondence between some operator O in the field 
theory and its dual bulk field <p in supergravity 

bulk field {supergravity) ~ operator {field theory). 

In particular, 

1. the dilaton field <j> is dual to the O = —£, where C 
is the Lagrangian density, 

2. the gauge field A a is dual to the i?-charge current 

3. the metric is dual to the stress-energy tensor in the 
sense that Z iD [g° u ] = e iScl ^^l 

More generally, the AdS/CFT correspondence states 
that 

Z AD [J] = J d(t>e lS+ S ldixJO = e lS "' 7 (9.3) 

and the partition function of the field theory Z[J], where 
J is the source coupled to the operator O. Differentiating 
Eq. (|9.3p with respect to J gives us various correlation 
functions. 

It is well known from the spectral theorem that the 
poles of the retarded Green's function of a wave equation 
coincide with the normal modes of the wave functions un- 
der appropriate boundary conditions. Thus, quasinormal 
modes of some gravitational background are naturally 
poles of the correlation functions. At the same time we 
can see that correlation functions in the field theory are 
connected through Eq. (|9.3[) with derivatives of the clas- 
sical action on the gravity side. Thus, it does not come 
as a surprise that the quasinormal modes of an asymp- 
totically AdS gravitational background coincide with the 
poles of the correlation functions 3 - 3 ^. When the field the- 
ory is at zero temperature, the dual gravitational back- 
ground is described by some regular asymptotically AdS 
metric. If the quantum field theory is at finite tempera- 
ture, the gravitational background must have a horizon, 
and the Hawking temperature of the horizon corresponds 
to the temperature in the dual field theory. A supposi- 
tion about this was done by Horowitz and Hubeny who 
suggested that the quasinormal modes of the large D- 
dimensional asymptotically AdS black holes are poles of 
the retarded Green's functions in the dual conformal field 
theory in {D — 1) dimensions 30 . As the calculations of 
the Green's functions directly on the CFT side in four 
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and higher dimensions are complicated, no proof of this 
QNM interpretation was found. Nevertheless, soon af- 
ter Birmingham, Sachs and Solodukhin 336 > showed that 
for the (2 + l)-dimensional asymptotically AdS (BTZ) 
black hole, the quasinormal modes, which are given by 
the relation^. 



uj = ±q — ATtTi{n + 1), 



(9.4) 



exactly coincide with the poles of the two-point retarded 
Green's function, 

.2 Jl 

G R = A- 



* 1 



i{u - q) 
4ttT 



4tt 2 



+ V 1- 



i(ui + q) 
4vrT 



(9.5) 



in two-dimensional CFT. This was the first proof of the 
AdS/CFT interpretation of QNMo 336 i 337 . 

Currently, the variety of gravitational backgrounds 
considered in string theory literature is quite large, be- 
cause of considerable interest in a number of condensed 
matter phenomena at strong coupling which can be mod- 
eled through the AdS/CFT correspondence. These are 
interesting problems because there are plenty of con- 
densed matter strong coupling phenomena (which in- 
cludes, for instance, strongly correlated electrons) which 
could be engineered in a laboratory. This area of re- 
search is now called the AdS/condensed matter the- 
ory (AdS/CMT) correspondence and one example of 
AdS/CMT problems, the holographic superconductor, 
we discuss later. Meanwhile, we mention here a recent 
review devoted to the AdS/CMT correspondence 3 - 3 -^ 

Next we discuss only the most universal properties of 
the AdS/CFT applications of QNMs which are valid for 
a wide class of dual gravitational backgrounds, thereby, 
which describe the hydrodynamic regime of strongly cou- 
pled conformal field theories. 



C. Universality of the hydrodynamic regime 

A great advantage that perturbations of black holes 
give through the AdS/CFT correspondence is the de- 
scription of the large distance collective behavior of the 
quark-gluon plasmas at strong coupling and, more gen- 
erally, the hydrodynamic regime of the dual field theory. 
In the regime of long- wavelength perturbations hydrody- 
namics can be described by the energy-momentum tensor 
of a system, 



T af3 = (e + P)i 



Pg 



a/3 



(9.6) 



where <r a P is proportional to derivatives of local temper- 
ature T and four- velocity u Q , is called the dissipative part 
of T a " . The equation of motion is the conservation law 



and, if there are other conserved currents, an extra con- 
servation law is added 



where 



f -o 



j a = pu a - V(g afi + u a uP)p. 



(9.8) 



(9.9) 



In the frame of reference in which the fluid is at rest, 
gives the law of diffusion 



this 



— > 

3 



-VVp. 



(9.10) 



Here, T> is the constant of diffusion. When consider- 
ing hydrodynamical processes related to the behavior of 
a quark-gluon plasma, one can distinguish three types of 
perturbations of the energy-momentum tensor which can 
be treated within independent sets of equations (in anal- 
ogy, for instance, with axial and polar perturbations of 
black holes). These modes are as follows 29 : 

• The shear mode, corresponding to perturbations 
of the T 01 , T 02 , T 31 , and T 32 components. This 
mode is responsible for translational fluctuations 
of the "cross section" of the system, a flow, also 
allowing, in the general case, for small rotations. 
In the frequency domain, the shear channel has the 
so-called hydrodynamic mode which is purely imag- 
inary and, in the linear approximation, is propor- 
tional to q>mm, 



-iq 



(9.11) 



where q is momentum normalized by the tempera- 
ture, q = k/2i:T. 

• The sound mode corresponding to perturbations of 
the T 00 , T 03 , and T 33 components. In simple sys- 
tems this mode is responsible for fluctuations in 
the direction of the flow, thereby generating fluc- 
tuations of energy density, i. e., sound waves. The 
sound mode usually has the following forr" 3 -^ 1 • 



u) = c s q - ijq 



(9.12) 



where c s is the speed of sound and 7 is some con- 
stant that depends on the energy density, pressure, 
shear, and bulk viscosities of the system 2 ^. 

• The bulk mode corresponds to perturbations of the 
T 12 component. Because of the conformality of the 
theory, the viscosity of the bulk mode equals zero. 
Therefore, we shall not discuss this mode in detail. 

Following 3 -! 2 . - ^!, we discuss two relatively easy and at 
the same time universal examples. 

Consider a quite general class of gravitational back- 
grounds given by the metric 



-■a/3 

,0 







(9.7) 



ds = g tt dt 



g rr dr 2 + g x 



(9.13) 
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The metric coefficients can be arbitrary, assuming oniy 
the following behavior near the event horizon r = r + : 



fjti 



-jt(r-r+), g r 



7r 



(9.14) 



Here jt & n d 7r are constants. As a particular example 
we consider the process of diffusion of the conserved R- 
charge in the J\f = 4 super Yang-Mills theory. The field 
equations for the gauge field dual to the conserved cur- 
rent have the form 



(9.15) 



where g e ff(r) is the effective gauge coupling. The con- 
served current j a , which is associated with the wave 
equation (|9.15[) . obeys the equation 



d a j a = 0. 



(9.16) 



Using the membrane paradigm, Kovtun, Son, and 
Starinets 342 showed that 



and found the diffusion coefficient to be 



(9.17) 



V = 



9xx9 e ffV~9tt9r 



(r+) / dr — (r). 



Indeed, the current j a can be written in the form 



(9.18) 



3 



(9.19) 



where is the radius of the stretched horizon in the 
membrane description. Then, it can be shown that in 
the radial gauged 3 - 

F ir = Vlr7ToF u (r -r+)-\ (9.20) 



F u « -diA t , (A r = 0). (9.21) 

The wave equation for the ^(-component of the gauge 
potential is 



8rW = gg rr 9 u 9rA t ) = 0. 



(9.22) 



The solution to the above equation (I9.22j) which vanishes 
at infinity has the form 



A t (r) = C 



°° d ,9tt(r')9rr(r') 



From the above equation (|9.23[) it follows that 



9tt9r 



-(r+) 



, 9tt9rr , v 

ar — Irj 



(9.23) 



(9.24) 



Then using the expressions for j and j 1 , one can find 
Fick's law (I9.10|) . From the above it follows that on the 
gravity side of the duality, independent of the details of 



the background metric, there is always the hydrodynamic 
quasinormal mode 



lu = -iVk 1 



(9.25) 



Another example concerns the calculation of the ratio 
of the shear viscosity to the density of entropy Ty/s 3 ^ 3 -. 
Following Buchel and Liu 3 ^ 3 - and Son and Starinets^, 
we perform the following boost of the metric f|9. 13[) : 

r = r', (9.26) 



t' + vy' . . 
t = = = « t' + vy', 



V? 



y = 



vt', 



(9.27) 



(9.28) 



Xi = x\. (9.29) 

Then, the metric (|9.13[) transforms to the following form: 

p 

ds 2 = g tt dt 2 + g„.dr n + g xx (r) ^(cfe'*) 2 + 

i=l 

+2v{g tt + g xx )dt'dy'. (9.30) 

The above metric is a k = perturbation, so that 

At = vg xx {g tt + g xx ) 

can be considered as a gauge potential, similar to the 
one in (|9 . 23|) . The potential vanishes at infinity and 
Eq. (|9.24[) now has the form 



A, 



1 + 9 xx 9tt , 



|r->r+0 



Ftr ' + d r {g**g tt y 
Then, the coefficient of diffusion is 



9xx(r+) 
7o 



V 



1 



4vrT 



(9.31) 



(9.32) 



Using the thermodynamic relations under zero chemical 
potential 



V = 



V 



p 1 



e + P = Ts, 



one finds 3 ^ 



j] 1 

S 47T 



(9.33) 



(9.34) 



Thus, for all field theories at finite temperature, which 
can be described by some dual gravitational background, 
the ratio of the shear viscosity r] to the entropy density s 
is always h/Ankf, to leading order in the 't Hooft coupling, 
where fcb is the Boltzmann constant. The next correction 
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in the series of powers of the inverse 't Hooft coupling is 
positivo 345 ! 346 , 

; = i( 1 + w)' « 3 >^ < M5 » 

This suggested a conjecture that for all systems that 
can be obtained from relativistic quantum field theory un- 
der finite temperature and chemical potential, the ratio 
n/s cannot be less than There are now various 

counterexamples to the above conjecture, showing that 
f]/ s ^ 1/47T. For example, this happens when consider- 
ing higher derivative quantum correction a 347 ' 348 . 

We now see what value of the viscosity/entropy ratio 
would produce the weakly coupled theory when naively 
extrapolated to a strong coupling regime. Viscosity, be- 
ing the measure of the diffusion of the momentum, is 
evidently larger for a larger temperature T and is de- 
creasing when the coupling A is weaker (because the free 
way of particles is longer which makes an easier transfer 
of momentum). More accurately, for example, for the 
\<j) A theory, one has2£ 

?7~S,s~T 3 . (9.36) 

Thus, 

- ~ A~ 2 . (9.37) 
s 

For the strong coupling regime A ~ 1 and 77/s ~ 1, which 
is one order more than the AdS/CFT value 1/47T. Ex- 
perimental data suggests that n/s is small, which means 
that the quark-gluon plasma is in the strongly coupled 
regime. The data extracted from relatively recent ex- 
periments with collisions of heavy atomic nuclei in the 
Relativistic Heavy Ion Collider— i 3 ^ say that if one mod- 
els the quark-gluon plasma by the Navier-Stokes equation 
with parameters fixed to fit the experimental data, then 
the ratio n/s must be indeed about one order less than 
unity and is close to 1/Att— This gives optimism to 
string theorists because the observation of the univer- 
sal n/s ratio in the RHIC or the LHC might be an ex- 
perimental confirmation of the AdS / CFT correspondence 
and direct connection with a real experiment, which was 
always lacking in string theory 3 ^- 3 ^ 2 -. 

D. Holographic superconductor 

In this subsection we showed how the WKB formula 
(I3.38p . used for finding QNMs, can be quite unexpect- 
edly used for estimation of the conductivity of super- 
conductors constructed in the spirit of the AdS/CFT 
correspondenc e 363 ! 364 . 

The effect of superconductivity consists of vanishing 
of the electrical resistivity of some metals (superconduc- 
tors) below some critical temperature T c . In addition, the 
magnetic field is expelled from such a superconductor at 



low temperature (Meissner effect), so that the supercon- 
ductor is also a perfect diamagnetic. The phenomenolog- 
ical theory of superconductivity by Bardeen, Cooper and 
Schrieffcr (BCS) 365 suggested that pairs of electrons with 
opposite spin can form a bound state by interacting with 
phonons. This bound state is a boson, called the Cooper 
pair, and below the critical temperature these pairs con- 
dense, thereby inducing superconductivity. However, the 
BCS theory works well for superconductors in the weak 
coupling regime, i. e., for those superconductors whose 
critical temperature T c is quite low. Superconductors for 
which T c is high are expected to involve strongly corre- 
lated electrons. For a description of this class of super- 
conductors the holography could provide some insight. 

On the gravity side of the gauge/gravity duality the 
superconductivity can be modeled by an asymptotically 
AdS black hole with scalar hair: A nonzero condensate is 
"represented" by the black hole "hair" . The field theory 
temperature is again dual to the Hawking temperature 
of the black hole so that one needs a black hole that 
would possess scalar hair at low temperature and would 
loose it at T > T c . Gubser showed that a charged black 
hole with charged scalar field around it satisfies these 
rcquirement o 363 ' 366 . Thus, the Lagrangian has the form 

Z = R+-^-\F^F^-\\7^-iqA<p\ 2 -Vm, (9-38) 

where ip is the scalar field, F^ v is the strength tensor 
of the electromagnetic field, q is the scalar field's charge 
and A is the vector potential (F = dA). The cosmological 
constant is —3/L 2 . 

According to the gauge/gravity dictionary the conduc- 
tivity can be found by solving for fluctuations of the 
Maxwell field. Therefore, further we study fluctuations 
of the potential A x . 

The plane symmetric solution can be written in the 
following general form: 

ds 2 = -g{r)e- x{r) dt 2 + ~7~T + r 2 {dx 2 + dy 2 ), (9.39) 



A = 4>{r) dt, = ip(r), (9.40) 

(9.41) 

x/2 

Using a new radial variable dz = s -^-dr, at large r, dz s» 

dr/r 2 , and we choose the constant of integration so that 
z = —l/r. The horizon is located at z = —00. Then, 
Eq. (|9.41|) has the following wavelike form: 

-A XjZZ + V{z)A x =u 2 A x , (9.42) 

where the effective potential i a 364 i 367 

V(z)=# 2 r + 2,V^]. (9.43) 
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Now, an incoming wave from the right will be partly 
transmitted and partly reflected by the potential barrier. 
The transmitted wave is purely ingoing at the horizon 
and the reflected wave satisfies the scattering boundary 
conditions at z — > oo, and will obey, at the same time, 
the Dirichlet boundary condition at z = 0. Thus, the 
scattering boundary conditions for z > are 



where 



A x — e 



Tie 1 ' 



and at the event horizon 



A x = Te 



-oo, 



(9.44) 



(9.45) 



where 1Z and T are reflection and transmission coeffi- 
cients. Then, one has 

A x (0) = l + n, A XtZ (0) = -icj(l-U). (9.46) 

As shown ir*2£& 



A, 



A 



(i) 



(9.47) 



Within the gauge/gravity vocabulary the limit of the 
electric field and the first subleading term on the bound- 
ary are dual to the field strength and the induced current, 
respectively, 



E,. 



-A® 

^x ' 



" X — -TL^ 



Thus, we find the conductivity 



a{uj) 



i A 



(1) 



"A x oy 



Since A x 1] = -A^O), then, 
cr(w) 



l-IZ 
1 + 1Z 



(9.49) 



(9.50) 



The above boundary conditions (|9.44[) . (|9.45p are noth- 
ing but the standard scattering boundary conditions for 
finding the S'-matrix. The effective potential has the dis- 
tinctive form of the potential barrier, so that the WKB 
approach can be applied for finding TZ and We 
note that as the wave energy (or frequency) uj was real, 
the first order WKB values for TZ and T are real 3 - 7 ^ and 
T 2 + TZ 2 = 1. Next, one can distinguish two qualita- 
tively different cases: First, when uj 2 is much less than 
the maximum of the effective potential Vq uj 2 <C Vq, and 
second when w 2 ~ Vq. When uj 2 ^> Vq, usually the reflec- 
tion coefficient 7Z decreases too quickly with u>, so that 
a reaches its maximal value (unity) even at moderate 
uj > Vq. Therefore, the case ui 2 ~ Vq also works well 
for large uj and one can use the WKB formula (|3.38p . for 
scattering around black holes related in (see Sec. IIII.Dl in 
this review), which gives 



(1/2))) 



-1/2 



uo 2 ~ Vq, (9.51) 



Xuj 2 -Vq) 



■ A 2 + A 3 



(9.52) 



Further and extensive literature on this subject con- 
siders various generalizations of the above example. 
This mainly includes the addition of a magnetic field 
and various corrections, such as higher curvature 
correction !] 364 ' 371 " — . 



X. CONCLUSIONS 

This review should be considered as an introduc- 
tion into the analysis of quasinormal modes of various 
black holes which are studied in astrophysics, higher- 
dimensional gravity, and string theory. Currently three 
comprehensive reviews on quasinormal modes 7 - - — , discuss 
in detail observational aspects of quasinormal ringing. 
We concentrated on the issues which were not reviewed 
or only briefly touched on there: methods of calculations 
of QNMs, stability of black holes and branes, late-time 
tails, the holographic superconductor, etc. Therefore, we 
sketched rather than discussed in detail some questions, 
referring the reader to more specialized literature. 
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